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Abstract 

In this thesis we consider several aspects of general relativity relating to 
exact solutions of the Einstein equations. In the first part gravitational 
plane waves in the Rosen form are investigated, and we develop a formalism 
for writing down any arbitrary polarisation in this form. In addition to 
this we have extended this algorithm to an arbitrary number of dimensions, 
and have written down an explicit solution for a circularly polarized Rosen 
wave. In the second part a particular, ultra-local limit along an arbitrary 
timelike geodesic in any spacetime is constructed, in close analogy with the 
well-known lightlike Penrose limit. This limit results in a Bianchi type I 
spacetime. The properties of these spacetimes are examined in the context 
of this limit, including the Einstein equations, stress-energy conservation and 
Raychaudhuri equation. Furthermore the conditions for the Bianchi type I 
spacetime to be diagonal are explicitly set forward, and the effect of the limit 
on the matter content of a spacetime are examined. 
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Notation 



• Unless otherwise specified c = 1, G^ = 1- 

• We work with signature (—1, 1, 1, 1). 

• Throughout V^ is used to denote the four dimensional Laplacian, with 
•^^^V^ the three dimensional Laplacian. Likewise when clarification is 
needed the dimension will be denoted in the same manner, e.g. ^^^R. 

Throughout this thesis we consider tensors in a number of different dimen- 
sions. To help clarify which number of dimensions we are working with in 
general: 

• lower-case letters from the start of the alphabet {a, 6, c, d} will be used 
for indices over four dimensions. 

• lower-case letters from the middle of the alphabet {i,j, A;,/} will indi- 
cate indices over three dimensions. 

• upper-case letters from the start of the alphabet {A, B, C, D} will de- 
note indices over two dimensions. 

Further we will use {a, (3, 7, 6} for labelling of an orthonormal basis, and we 
leave it for the reader to determine the dimension based on the tensor indices 
used. 
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Preface 



This thesis focuses on two separate topics both of which are related to phys- 
ically interesting exact solutions in classical general relativity. 

The first research problem concerns polarisations of strong-field gravitational 
plane waves. This is one case where an exact solution is simple enough to 
be mathematically tractable but, although highly idealised, may still capture 
aspects that are of physical relevance. 

As with many well-studied metrics in general relativity, more than one 
way has been discovered to represent these gravitational waves, known re- 
spectively as the Rosen and the Brinkmann forms. These two forms have 
their respective advantages, depending on the scenario considered, as would 
be expected. However it is important to know, in theory, how to construct 
solutions for each form. When there are two ways of representing the same 
phenomena, and it appears much harder, or impossible, to formulate solu- 
tions in one form than the other, then it may indicate that work in this 
respect is incomplete. 

In the case of gravitational waves, it was trivial to represent arbitrary 
polarisations in the Brinkmann form, while it was unclear how to construct 
arbitrary polarisations in the Rosen form. As such we investigated: 

• Is there a clear way to construct arbitrary polarisations in the Rosen 
form? In particular is there a simple algorithmic construction to do 
this? 

• Is it possible to write down particularly useful polarisations, such as 
circular polarisation in the Rosen form? 

These questions have been answered positively, and a general method for 
constructing such polarisation modes has been developed. 
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IV 



One context in which plane waves, and their generahsations, pp waves, are 
studied is as the result of a Penrose limit. This well-known limit approxi- 
mates the spacetime along any arbitrary null geodesic in any spacetime. We 
attempted to find a similar, timelike, limit. As such the questions investi- 
gated in this thesis are: 

• Is it possible to develop a similar limit for timelike geodesies, and what 
is the result of this limit? 

• In what way does this limit share the properties of the Penrose limit? 

• What happens to any matter content when taking the limit? 

We have developed a limit that appears to be a close timelike analogue 
to the Penrose limit, although our construction does not share all of the 
important properties of Penrose limits. It is also found to be a type of ultra- 
local limit, whereby the evolution of neighbouring points decouples. 

The result of this limit was a Bianchi type I spacetime. This is another 
well-studied exact solution which is simple enough to be mathematically 
tractable. It has most frequently been studied as a basic cosmological model. 
This has led to most studies making context-appropriate assumptions that in 
the context of this thesis should not be made. In particular most studies have 
worked with a diagonal metric, and mostly examined vacuum, perfect fiuid 
or dust scenarios. In counterpoint, in our context we found the most general 
scenario for the stress-energy was that the limit would erase information 
about the fiux, but without more stringent conditions for general spacetimes. 
Therefore some further questions were raised: 

• What can we say about Bianchi type I spacetimes whilst keeping full 
generality? 

• What can we say about when it is appropriate to have a diagonal 
metric? 

We have provided a set of necessary and sufficient conditions for the met- 
ric to be diagonal. We have furthermore made a number of calculations 
relating to the Einstein equations, covariant conservation and the Raychaud- 
huri equation in general Bianchi type I spacetimes, simplifying as much as 
possible, and interpreting the results. 



Arrangement of this thesis 



This thesis has been arranged into four main chapters plus conclusions. The 
first is a general introduction to some aspects of general relativity, many 
of which are used later in the research sections of the thesis. Much of the 
information in this chapter could be found in a good introductory textbook 
in general relativity. 

The second chapter introduces and reviews the topic of strong-field gravi- 
tational waves before going on to examine the problem of constructing strong- 
field polarisations of the Rosen form. From section ( 12. 3 p onwards the infor- 
mation presented is primarily original research, results from which have been 
presented in [T8| [T9] . 

The third chapter forms another background section providing informa- 
tion on some topics of particular relevance to chapter four, which develops 
and explores a Penrose-like limit. Unlike the first chapter, much of this infor- 
mation (with the exception of the discussion of the FLRW metric) would not 
be included in a basic textbook, but is more advanced information familiar 
to many researchers working in the relevant areas of general relativity. 

Chapter four consists of the construction of, and investigation into, an 
ultra- local, Penrose-like limit. This chapter primarily consists of original 
work, some results of which have been published in |2U| . 

The appendix provides the abstracts and publication details for both 
journal papers and the conference article that have resulted from the research 
presented in this thesis, all of which are coauthored with Prof. Matt Visser. 
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r CHAPTER 1 

Introduction to general relativity 



Here we will very briefly describe some of the important aspects of general 
relativity, focusing on techniques and results that will be important to the 
research presented in this thesis. 

General relativity is a theory of gravity first formulated in around 1915, 
primarily by Einstein, with significant contributions from Hilbert, and which 
has been developed extensively in the intervening decades. In general rel- 
ativity the force of gravity is due to spacetime curving in the presence of 
matter. The natural mathematical language for considering curved spaces is 
that of tensors and differential geometry. 

1.1. Geometrical tensors 

Spacetime is described by the four-dimensional metric tensor, Qat, with the 
associated invariant interval 

ds^ = gabdx''dx''. (1.1) 

The metric tensor describes the infinitesimal notion of distance in various 
directions. The important geometrical quantities in general relativity are 
formed from the metric tensor and its derivatives. 

In particular, when considering the acceleration due to gravity it is natural 
to consider tensors formed from up to second derivatives of the metric. There 
are a limited number of these. 

Firstly, it is helpful to introduce the connexion/Christoffel symbol {not a 
tensor), given by 

r%c = - g""^ {gdh,c + gdc,b - gtcd) ■ (1-2) 

The Christoffel symbol acts to propagate a vector infinitesimally from one 
point to another, and can be defined in terms of the covariant derivative of 
a tensor 

Vf,T'^ = 9feT'^ + r^5T^ (1.3) 

Now the Riemann tensor can be formed from the Christoffel symbols and 
their derivatives 

R"'bcd = dc T'^bd — dd Vbc + r"ec ^'^bd — ^"'ed ^^bc- (1-4) 



The Riemann tensor describes the curvature of the spacetime, and can be 
associated with tidal forces on nearby test particles. 

The Riemann tensor can be formed from two simpler tensors; 

Rabcd = Cabcd + " {dac Rbd + dbd Rac — dad Rbc — Qbc Rad) 

- g^ {gacQbd - gadQbc) , (1-5) 

where Cabcd is known as the Weyl tensor, defined by the above equation, the 
significance of which is that it is the curvature which remains present in a 
vacuum, and 

Rab = R'^acb] (1-6) 

R = g'''Rab, (1.7) 

are the Ricci tensor and Ricci scalar, respectively. 

Finally, the Einstein tensor is constructed from the Ricci tensor and Ricci 
scalar. 

Gab = Rab - -^Rgab ■ (1-8) 

1.2. Geodesics 

Geodesies are the paths, often denoted by 7, followed by test particles, influ- 
enced only by the curve of spacetime. Here by a test particle we mean the 
mass of the particle does not affect the spacetime through which it travels. 
There are three sorts: 

• spacelike, followed by hypothetical faster-than-light particles (tachyons) 
and therefore of lesser physical interest 

• lightlike, also called null, followed by particles travelling at the speed 
of light (luxons) 

• timelike, followed by slower than light particles (tardyons/bradyons) 

Understanding geodesics is of prime importance when trying to understand 
the physics of general relativity. 

In familiar Newtonian physics, the path of a free particle is described by 

^^0, (1.9) 



indicating no forces act on the particle. Geodesic motion is similar; no forces 
act on the particle and it simply falls under the influence of the spacetime 
metric. In Newtonian mechanics, a free particle follows a straight line; like- 
wise geodesies are the straightest possible paths, which are also the paths of 
shortest length in the appropriate sense. 

To describe this straightest path we take a tangent vector at a point, T", 
and require that the tangent vector can be parallel transported so that the 
transported vector is proportional to the original, 

T^VaT^ (xTK (1.10) 

Consider an arbitrary parameterisation of a curve x°(A), which has the 

tangent -tt-. Insert these into equation fll.lOp to give the geodesic equation 
dA 

with /(A) some arbitrary function. It is always possible to write this in 

an affine parameterisation by an appropriate choice of A, which simplifies 

equation (11. lip to 

^2 a ^ b ^ c 

Lightlike geodesies have, for any parameterisation 

dx'^ dx^ 

For timelike geodesies, an affine parameter can be chosen to satisfy 

dx" dx* _, /, , .N 



1.3. Coordinate independence 

For every scenario there are many possible representations, many possible 
choices of coordinates. The choice of coordinates selected might have some 
advantage in ease of representing a particular feature of the spacetime. 



As a simple example of this, consider the Schwarzschild solution, a non- 
rotating black hole. The standard choice of coordinates, often called Schwarz- 
schild coordinates or curvature coordinates is 

ds^ = - (l - ^) dt^ + -^^ dr^ + r^ {d9' + sin2(^) d^^) . (I.15) 

However, numerous other coordinate choices have been used over the years, 
such as for instance the Painleve-Gullstrand coordinates 

ds^ = - M —J dP+2J^dfdi+df^+f^{de^ + sm^{e)d(j)'^). (1.16) 

Two metrics represent the same spacetime if there exists a transformation 
between the two. Often explicitly finding this transform is the best way to 
show two metrics are equivalent. 

There is a transformation of the time coordinate between these two forms. 
The major physical difference between these forms is that at the event horizon 
at r = 2M, in Schwarzschild coordinates the dr component becomes infinite, 
resulting in what is known as a coordinate singularity, whereas the Painleve- 
Gullstrand coordinates can be extended across the horizon. However the 
Schwarzschild form has a simpler diagonal form, and the meaning of the 
time coordinate in the two forms is different. 

Since the best coordinates to use depends on what properties one wants 
to investigate, and in many cases different choices will be best for different 
properties, it good to understand (at least in principle) the basic construction 
of important features in different commonly used systems. 



Often coordinates are chosen because they represent what is seen by an 
important class of observers. For instance, consider synchronous coordinates, 
which are of the form 

ds2 = _dt2 + g..(^t, x) da;Mx^'. (1.17) 

This means that paths with fixed spatial coordinates, x* = x^o, and varying t 
are geodesies. The coordinate t represents the proper time of these observers. 
As such synchronous coordinates describe what is seen by freely falling ob- 
servers. Such coordinates often cannot be globally extended as there are 
coordinate singularities where two such paths cross. Thus it is not always 
possible to pick such a system globally, but in theory is always possible to 
do locally (though technically difficult in practice). 



1.4. Tetrads 



In any coordinate system it is possible to pick what is known as a tetrad 
frame (also referred to as a vierbein, orthonormal basis, or non-coordinate 
basis). 

To construct such a basis, at any point pick four linear independent vec- 
tors, e^Q, (here a is merely a way to number which vector we are talking 
about) that are orthonormal in the sense 



9ah^ aG 13 — Va/S, 



where 



Val3 



-10 
10 
10 
1 



:i.i8) 



;i.i9) 



is the metric of Minkowski space, appropriate for special relativity. 

The matrix formed from these contravariant vectors has an inverse matrix 
which is formed from the covariant vectors e^a- Given this set up, tetrads 
form an orthonormal basis to the tangent space at a point. This implies 



dab 



3/3 



bVa/S- 



;i.2o) 



Such objects are also useful in other dimensions, where in general these 
are known as vielbeins. Particularly useful are three dimensions (triads or 
dreibeins) and two dimensions (zweibeins). 

Tetrads are useful for comparing coordinate-based calculations to what 
would be observed in a given scenario, as each observer essentially carries 
a tetrad basis (their standard notion of space and time) with them. In 
particular to make useful comparisons to the matter /energy that would be 
observed, one should use a tetrad basis. Objects calculated in a tetrad frame 
(such as the stress-energy) will be denoted by hats over the indices. 

1.5. Matter 

In addition to tensorial descriptions of the spacetime and its properties we 
want tensorial descriptions of any matter that may be present in the space- 
time. This is achieved with the stress-energy tensor: 



P 



Z 



ab 



1.211 



TTi 



with p the density, Fi the flux of energy, and tTjj the stress tensor, which 
describes the pressures and stresses present. 

Of course, the form of the stress-energy is coordinate dependent. For 
instance, an observer passing through a field of dust will notice a flux, whereas 
one moving with the dust will not. Several useful types of matter to examine 
are: 



1. Vacuum 

Vacuum has T"-^ = 0, independent of the coordinates used. 

2. Dust 

Dust is a from of matter that is without pressure, given by 

with V"" timelike. 

In a comoving orthonormal frame it is simply given by 



;i.22) 



T 



ab 



p 









;i.23) 



3. Null dust 



Null dust is the lightlike version of pressure-less matter, suitable for de- 
scribing photons and hypothetical speed-of-light particles. The stress- 
energy tensor can be decomposed into the product of a null vector with 
itself, similarly to the dust case. 



If we employ null coordinates 

u = z — t; V = z + t, 

there is a frame where the only component is T^u- 
4. Perfect fluid 



;i.24) 



;i.25) 



A perfect fluid has density and pressure, but no viscosity, and hence 
cannot support a shear, and additionally all pressures are the same. A 
perfect fluid is described by 



T'''' = {p + p)V''V^ + pg''\ 
so in a comoving orthonormal frame it takes the form 



db 



p 








" 





p 














p 














p _ 



T 



5. Scalar field 

The stress-energy of a scalar field can be expressed as 

1 



T,, = 0,a0,,--(7„, {(V0)=^ + \/(0)} 



;i.26) 



;i.27) 



:i.28) 



Examples of (hypothetical) fundamental scalar fields are the Higgs field 
and dilatons. 

6. Electromagnetic field 

The electromagnetic stress-energy tensor is given by 



T, 



ah 



\ {E^ + B^) Ex B 

Ex B EiEj + Bi Bj -^{E^ + B^)5ij 



;i.29) 



containing the familiar Maxwell stress tensor, and the Poynting vector 
as flux. 

It can sometimes be convenient to express fields using the vector po- 
tential, A"-, 

A''{t,x) = ((l){t,x); A{t,x)) , (1.30) 



:i.3ii 



which describes the electric and magnetic fields through 

dA 



B = V X A; E = -V0 - 



dt 



We can express the electromagnetic field tensor as 

-Tab = ^b;a ~ ^a;b} 

or explicitly, 



;i.32) 



ab 



Ex Ey Ez 

—Ex Bz —By 

—Ey —Bz Bx 

— Ey B,, — B^. 



(1.33) 
In terms of the electromagnetic field tensor, the stress-energy is given 

by 

Tab = Fac g'^Edb - \ gab {Fef F'f) . (1.34) 

Note the trace of the stress-energy vanishes, g"" Tab = 0. 

1.6. Einstein equations 

This brings us to the central equations of general relativity, the Einstein 
equations: 

Gab = 8nGNTab. (1.35) 

These equations link the presence of matter to the curvature present in the 
spacetime, and determine how that curvature will evolve. 

Although simple in appearance these are ten second-order non-linear par- 
tial differential equations in the general case, which have only been solved 
exactly in situations of high symmetry [32| [62] . 

However, these ten equations are not independent, and are constrained 
by the Bianchi identity, which, in the context of general relativity (rather 
than abstract differential geometry in general), is 



VaG- 



ab 



0. 



;i.36) 



This constraint means that only six of the ten of the Einstein equations are 
independent, and we thus have coordinate freedom to reduce the equations 
down to six. 



The Bianchi identity can also be regarded as a conservation relation on 
the stress-energy tensor, 

VaT"'' = 0, (1.37) 

as a consequence of the Einstein equations. 

For vacuum the Einstein equations reduce to a statement about the Ricci 
tensor as 

Rab--R9ab = => Rab = 0, (1.38) 

which gives the vacuum equations. 

1.7. Weak-field gravity 

Although solving the general equations in situations without significant sym- 
metry is difficult, it is possible to explore the weak-field, linearised results 
more generally. This was one aspect that was investigated soon after the 
development of general relativity, primarily to predict differences from New- 
tonian gravity. 

When a gravitational ffeld is weak it is expected that the spacetime metric 
should be close to the description of special relativity; 

ds^ = r]abdx''dx\ (1.39) 

So the linearised theory of of gravity is obtained by using the metric: 

gab = flab + hab\ kab <t: I, (1.40) 

and neglecting all second order or higher terms when calculating relevant 
physical variables, for instance 

r^c = \ V"' {hdb,c + hdcb - hc,d) + 0{e), (1.41) 

Rab = \ {h\,bc + h\ac - V hab - kab) + 0(/l=^). (1.42) 

In addition we have gauge freedom to set 

h'b-^h6^\ =0. (1.43) 
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This is known variously as the Einstein/Hilbert/de Donder/Fock gauge, and 
siniphfies equation (ll.42p down to 

Rab = -lv'hab + 0{e), (1.44) 



2 



and consequently, 



Gab = -\ V [hab - \ hi]A + 0{y?). (1.45) 

So we have reduced the Einstein equations to simpler, linear differential equa- 
tions. 

Note the strongest gravitational fields in our Solar System are well within 
the regime where these approximations are a very good description of reality. 
This means that weak-field gravity has been very important for measurable 
predictions and applications. These include time dilation due to gravity wells, 
which is important in GPS applications and the precession of perihelion of 
Mercury that had been noted prior to the development of general relativity. 
Another important prediction is that gravity propagates through spacetime 
as a wave, direct evidence for which is still being sought. We will examine 
these weak-field waves in some detail. 

1.7.1. Weak-field gravitational waves 

From equation (ll.44p it is clear that the vacuum solutions for weak-field 
gravity are 

V^hab = 0. (1.46) 

This means that the metric obeys the wave equation in a vacuum, and the 
effects of changes in gravitational fields propagate via gravitational waves. A 
way to note the physical effect of these waves is to look at the (weak-field) 
Riemann tensor, 

Rabcd = ^ {had,bc + hij^^ad — hac,bd — hhd,ac} , (1-47) 



and so 



V^ Rabcd = 0. (1.48) 



We can investigate the properties of gravitational waves through their effects 
on test particles due to tidal forces apparent by the non-zero Riemann tensor. 
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Consider a plane wave travelling in the z direction, which we could write as 

ds^ = -dt^ + dz^ + (5^^ + hABit, z)) Ax^ dx^, (1.49) 

or alternatively, 

ds^ = dudw + ((5AB + /iAB(«))da;^da;-^, (1.50) 

where we are employing null coordinates: 

u = z — t; v = z + t. (1-51) 

Note there are several natural definitions for null coordinates that are com- 
monly seen, such as 

u = -^{z-t); v = -^{z + t), (1.52) 



72*^-''' 


-^(^- + *) 


u = t — z; 


V = t + z. 



(1.53) 

These minor differences lead to frequently seen differences in factors of ±1 
or ±2. 

For the wave given in equation (ll.50p the components of the Riemann 
tensor are Ruyuyi Ruxuxi and Ruxuy, so there is no motion of test particles in 
the z direction. Therefore we will consider a circle of test particles perpendic- 
ular to the 2-axis. It is found that the circle is distorted whilst maintaining 
the same area. 

To concisely describe the way that such a circle of particles is distorted 
we introduce polarisation modes. For gravitational waves the two modes are 
45° out of phase. They are referred to as the + and x polarisations. The 
effect of these are shown in figure (II. ip . Pure + polarisation can be written 
as 

ds^ = dudv + dx^ + dy^ + - f{u) [dx^ - dy"^) , (1.54) 

and the x polarisation as 

ds^ = dudv + dx^ + dy"^ + g{u) dxdy. (1.55) 

As superposition is possible due to the linearity of the wave equation, 
these two modes can simply be combined to create any weak-field gravita- 
tional wave, with some appropriate choice of f{u) and g{u): 

ds^ = dudv + dx^ + dy^ + ^-^ (dx^ - dy"^) + g{u) dxdy, (1.56) 
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Figure 1.1: The effect of the + (upper) and x (lower) polarized gravitational 
waves, on a circle of particles with the plane of the page perpendicular to the 
direction of propagation 



where in all cases f{u) ^ l,g{u) <^ 1. 

For instance, a circularly polarized wave can be constructed by 



ds^ = dudv + dx^ + dy^ + S I cos{ku] 

with 5* ^ 1 some constant. 

1.8. The Raychaudhuri equation 



(dx^ — dy"^ 



m{ku) dxdy j , (1.57) 



sm 



In general relativity the Raychaudhuri equation [Ml 165] is used to describe 
the motion of nearby bits of matter and whether they converge or not. There 
are two versions, for lightlike and timelike particles. For the timelike case 
one starts with a collection of nearby timelike vectors, u"" (a congruence). 
An example of this would be the 4-velocity of a fluid. Define the projection 
operator 



hab = gab + Ua Ub- 

Define the expansion tensor and scalar. 



:i.58) 



;i.59) 



13 



e = g''''eab = Vau\ (1.60) 

The expansion tensor measures the rate at which the volume of a small ball 
of fluid changes. The shear tensor and scalar are: 

(7ab = dab- ^hab9, (1.61) 

a' = ^ ^afe cr'^' > 0. (1.62) 

These measure the tendency for a small ball of fluid to be distorted into an 
ellipsoid shape. Now consider the vorticity tensor and scalar, 



^ab 



hacVK'^^hdb, (1.63) 



CO^ = ^UJabUJ^'>0. (1.64) 

These measure how the nearby vectors twist around each other. 

Note there exist some minor differences in commonly used notation. We 
are following the conventions of Hawking and Ellis [34j. 

The expansion, shear and vorticity tensors can also be viewed as a decom- 
position of VbUa into a symmetric traceless tensor (shear), symmetric tensor 
with trace (expansion), and an antisymmetric tensor (vorticity). That is, 

VbUa = u;ab + crab + -(^hab- -r^Ub. (1.65) 

3 ds 

Raychaudhuri's equation is 

^ = -R.b u^ u' + 2a;2 -2a'-U' + V. {^\ . (1.66) 

The lightlike case is similar, except that the definition of the projection 
operator hab as in equation f ll.SSp does not work in quite the same way, due 
to technical issues with constructing orthogonal vectors to a null ray. This 
is because the null ray has zero length, hence is orthogonal to itself. 

Due to the physically useful picture these terms provide it is often of 
interest to describe the form of the vorticity, shear and expansion tensors in 
particular spacetimes, especially if the Raychaudhuri equation simplifies in 
some way. The major use for this equation has been in the development of the 
singularity theorems [31]. It has also been of use in such diverse contexts as an 
astrophysical understanding of when cracks form in dense objects (see section 
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IV of [40] and the references therein), and of importance in understanding 
the gravitational analogues of the laws of thermodynamics. The lightlike case 
is particularly useful in examining gravitational lensing. Various extensions 
have been made to this equation in cases such as speculative theories of 
gravitation involving torsion. See f27\ 140] for reviews. 



1.9. Splitting of spacetime into space and time 

The Arnowitt, Deser and Misner (ADM) decomposition Pj STf is a way of 
splitting a general metric into propagating slices or "leaves". This splitting is 
referred to as a foliation. The spacetime metric is decomposed into the lapse 
function, N, shift vector, N\ and spatial metric Qij] 



gab 



Then the inverse metric is 



" NkN'' - Ar2 


Nk 1 


[ Nk 


9ij 



;i.67) 



^ab 



1 

Ar2 


Ar2 


iV2 


g'i -N'N^/N^ 



(1.68) 



This decomposition has been used in various attempts to construct a 
quantum theory of gravity. The construction is also useful for understanding 
one of the common concepts of mass in general relativity, the ADM mass. 

1.9.1. Extrinsic and intrinsic curvature 

There are two sorts of curvature that may be of mathematical and physical 
importance. General relativity primarily concerns itself with the intrinsic 
curvature (of four- dimensional spacetime), which is the curvature which can 
be measured by observers stuck inside the space itself, and is measured by 
the Riemann tensor. Another important way of looking at curvature is the 
extrinsic curvature, which is the curvature measured with respect to an em- 
bedding in a higher-dimensional space. The most notable simple difference 
between these two concepts is for an arbitrary ID space, which has zero in- 
trinsic curvature (the Riemann tensor is always zero in ID) but which can 
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be defined as curving by placing it in a 2D space and ascribing a radius of 
curvature to tlie various points on tlie curve. 

Let us make tliese two concepts precise: take a foliation of spacetime, and 
take a normal vector, Af, along some path, s, in a given slice. This normal 
projects a slice infinitesimally into the future slice. The extrinsic curvature, 
K, measures how the normal changes direction from one point on a slice to 
another. 

A' = ^. (1.69) 

as 

Although the intrinsic curvature is of most interest in general relativity, 
we saw above that it is possible to split 4D spacetime into 3D spatial slices 
propagating in a fourth, timelike, dimension. This means we can discuss 
the extrinsic curvature of the spatial slices with respect to the timelike co- 
ordinate. In fact, it is always possible to express intrinsic four dimensional 
curvature in terms of the extrinsic curvature and the intrinsic curvature of 
the spatial slices |47| . Depending on the particular spacetime considered this 
may be a very useful property. 

This can be particularly useful when one is working in the ADM view- 
point, as we may calculate the extrinsic curvature as 

K^, = ^ (iV.. + N,, - ^) , (1.70) 

where the covariant differentiation is with respect to the three dimensional 
metric, gij. 

1.9.2. Lagrangian formulation 

The Lagrangian for (vacuum) general relativity is 

£ = ^/-Wg W/?. (1.71) 

In the ADM viewpoint there is a simple expression for the Lagrangian 

C = (3)^1/2 ^ ^K,, K'^ -K^ + (3)i?) . (L72) 

This means that we have the simple interpretation that 

KijK'^-K^ (1.73) 

can be regarded as kinetic energy, and ^^^R as the potential energy. 
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1.10. Discussion 

In this chapter we have outhned the basics and some key important results 
from the theory of general relativity. We have been highly selective, focusing 
only on items that are relevant to the rest of this thesis. 

The Einstein equations link the evolution of spacetime with the matter 
content, and we can describe these spacetimes with a metric tensor. Further- 
more even in the case of weak-field gravity, general relativity makes some 
wildly different predictions form Newtonian gravity, such as the existence of 
gravitational waves. 

In addition we saw a useful way to split spacetime and discovered some 
nice ways to examine the behaviour of matter, through the geodesic equations 
and the Raychaudhuri equations. 



r CHAPTER 2 

Polarisations of the Rosen form 



In the introduction we looked at weak-field gravitational waves. Now we 
will consider their strong-field analogues, in the two common forms they are 
presented in, the Rosen and Brinkmann forms. We find in one of these forms 
it is a simple matter to write down any arbitrary polarisation, while in the 
other form it is not clear how to construct such polarisations. 

To address this issue we develop an algorithm to write down arbitrary po- 
larisations, filling this gap in our understanding of strong-field gravitational 
waves. The research in this chapter has been presented in [TSl [TU] . 

2.1. Strong-field gravitational waves 

We have a complete description of weak-field gravitational waves, which will 
include every gravitational wave we can potentially directly detect on Earth. 
However strong-field, non-linear, gravity is important in many astronomical 
events, which will produce gravitational radiation, and thus it is important 
to understand the strong-field gravitational waves. In general this is a very 
complex problem, so one aspect that is commonly investigated is the case 
with planar symmetry. While this is an unrealistic assumption, it can be 
hoped that some general features of strong-field gravitational waves may be 
understood by understanding this simplest case, and the planar case already 
demonstrates some interesting features and raises some important questions. 

Note that strong-field gravitational waves do not obey the wave equation: 
this is not a problem, and is not unique to gravitation, instead being a re- 
fiection on the assumptions which are necessary to obtain the wave equation. 
This means strong-field waves do not obey superposition, as the Einstein 
equations are not linear, and the investigation of colliding gravitational waves 
has been of considerable interest within the relativity community [T5l [32] . 

Progress was made in finding the strong-field analogues of weak-field grav- 
itational waves, in situations of high symmetry, by Einstein and Rosen [25j 
and independently by Brinkmann [12], in two different forms. We shall ex- 
amine both these forms of strong-field gravitational plane waves, starting by 
introducing a specific ansatz based on our knowledge of weak-field waves. 
Finally we shall show that the two forms are related by a coordinate trans- 
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form, and hence describe the same phenomenon. The respective advantages 
and disadvantages of the two forms will be mentioned. 

2.1.1. Rosen form 

We have found we can write weak-field waves as 

ds^ = dudv + {5ab + hAB{u))dx^ dx^ . (2.1) 

One ansatz for a strong field wave is simply to remove the restriction that 
^ab(^) is small, so 

ds^ = dudv + gAB{u)dx^ dx^ . (2.2) 

It is then a standard result [62j that the only non-zero component of the 
Ricci tensor is: 

Ruu = ~ j 2 ^ ^^^ ~ 4 ^ ^'^'^ ^ ^''^^ I ^ ^""' ^'^''^^ 

as i? = 0. 

Vacuum solutions for this metric are strong field gravitational waves in 
the Rosen form. Note that for weak-field waves hAsiu) was a completely 
arbitrary symmetric matrix function of u, while in the strong- field case we 
require Ruu = 0, that is, 

I /^ 9'ab - \ /^ 9'bc f"" 9'nA = 0, (2.4) 

so we have a restriction on qab- 

It is easily found, up to the usual symmetries, that for a wave travelling 
in the 2;-direction, the only non-zero components of the Riemann tensor are 
Ruxux, Ruxuy and Ruyuy, givcu by 1 18] 



RuAuB = ~ ] 2 ^^^ ~" 4 ^^'^ ^ ^^^ I ■ '^^■^'' 

Considering some of the properties of this form we note that it is explicitly 
translation invariant in the directions perpendicular to motion, as expected 
from the symmetries for a plane wave. However this form is unfortunately 
prone to coordinate singularities which indicate the formation of caustics of 
geodesies. 
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2.1.2. Brinkmann form 



Consider the metric 

ds^ = dudv + {Hab{u) x^ x^} dw^ + dx^ + dyl (2.6) 

This is the form for a plane wave in the Brinkmann form |[12j. The only 
nonzero component of the Ricci and Einstein tensors is given by 

Ruu = -- {Hxx{u) + Hyy{u)) = Guu, (2.7) 

again, with R = 0. 

Furthermore (up to the usual symmetries) the non-zero components of 
the Riemann tensor are 



1 

2 



J^uyuy „ ^yyV^)i l^-^j 



Ruxuy = —-^Hj,y{u). (2-10) 

Note that here, and with equation (12. 7p . the subscripts label components, 
and do not indicate derivatives. 

We see that this is the same pattern of non-zero components as for the 
Rosen form. This means that some of our intuition, developed with weak- 
field Rosen form, can be carried over to strong- field Brinkmann form |32] . 

Note this form is explicitly not translation invariant in the x and y direc- 
tions, unlike the Rosen form. However, this form has a distinct advantage in 
that it is not plagued by coordinate singularities in the same way that the 
Rosen form is. Thus we see that one form of strong-field gravitational waves 
better represents the symmetry of the system at the cost of the coordinates 
being inextendable beyond a certain point. 

Gravitational plane waves are then given by solving the algebraic vacuum 
equations 

H,,{U) = -Hyy{u), (2.11) 

or more explicitly, 

' Hju) HJu) 



Hab 



' X I 



HJu) -HJu) 



■ X 



(2.12) 



with HJu) and HJu) completely arbitrary (the reason for this choice of 
labelling for these functions will become clear shortly). 
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2.1.3. Transformation between the two forms 

To show these two forms really represent the same phenomenon, we will now 
demonstrate the coordinate transformation that takes the Brinkmann form 
into the Rosen form (see [SI E2])- 
Starting with the Rosen form 

ds"^ = dU dV + gAB{U) dy"^ dy^ , (2.13) 

we will transform to the Brinkmann form: 

ds^ = dudv + H^p{u) x°x^ du^ + dx^. (2.14) 

Begin by transforming the transverse coordinates as 

a;" = E°By'', (2.15) 

where we use x'^ for transverse coordinates in Brinkmann form, and y"^ for 
those in Rosen form, and E"" b obeys 

gAB{U)=E''AE^B^a^. (2.16) 

So it is a zweibein. Now 

dy^ = E^^ dx" + E^p x^dU, (2.17) 

and 

QAB d/ dy^ = E'^A E^B 5^p{E^^ dx^ + E^^ x^dU) {E^s dx^ + E^s x^dU). 

(2.18) 
Now simplifying 

QAB d/ dy^ = 6ap (dx" + E° A E^^ x^ dU) {dx^ + E^ b E^s x^ dU) , (2.19) 

and expanding out 

QAB d?/^ dy^ = 6ap (dx^dx'^ + E^b E^s x^ dUdx" 

+ E'^A E^^ x^ dUdx^ + E'^A E^^ x^ E^B E^s x^ dU^) . (2.20) 



This gives us a dU"^ term, and the desired form for the dx°' dx^ components, 
but with undesired dx'^dU terms. To complete the transformation we trans- 
form v. 

V = v + E^A E^s x" x^ = IgAB / y""; (2.21) 
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dV = dv-]^ QAB V^V"" -\ 9AB dy^ y^ - ^ QAB V^ dy^. (2.22) 

Now using equation ( 12.1 7p to write y^, dy^ in terms of x" and E'^a we see 
the dxdu terms cancel as well as producing new du"^ terms, giving us the 
desired Brinkmann form; 

ds^ = dUdV + df^ + Ec^i E'p x" x^ df/^. (2.23) 

So the desired transformation from one coordinate system to another is: 



u 


= u; 


1^ 


_ T^A a. 


V 


= v + E^aE^^x^x'' 



(2.24) 

It is thus seen that these two forms are related by a coordinate transform, 
hence represent the same phenomenon. Note that this transformation is 
singular at any point the Rosen form has a coordinate singularity. 

2.2. The pp waves 

There is a generalisation of the concept of gravitational plane waves; the pp 
waves. Consider the spacetime geometry [^ [SU EH E2] 

ds^ = dudv + H{u, X, y) du^ + dx^ + dy"^. (2.25) 

This is the Brinkmann form for a general pp spacetime. 

It is then a standard result that the only nonzero component of the Ricci 
tensor is 

Ruu = ~\{dlH{u,x,y) + dlH{u,x,y)] = Gun. (2.26) 

Where again i? = 0. Furthermore (up to the usual index symmetries) the 
only non-zero components of the Riemann tensor are of the form RuAuB- 
Specifically: 

Ruxux = --dlH{u,x,y)- (2.27) 

Ruxuy = -- 5x dy H{u, X, y); (2.28) 

Ruyuy = -ld'yH{u,x,y). (2.29) 
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This metric describes pp wave (plane-fronted waves with parallel propaga- 
tion) spacetimes, of which plane gravitational waves are a subset pi]. This 
spacetinie describes gravitational waves where the source is either electro- 
magnetic waves, waves of any hypothetical zero- mass particle, or null dust. 

One notable aspect of pp waves is that [8j all the scalar invariants vanish 
(such a spacetime is referred to as a vanishing scalar invariant (VSI) space- 
time). This means a pp spacetime cannot have a scalar curvature singularity 
(where the curvature invariants become infinite), though it may have other 
sorts of singularities. 

In addition to plane waves, some important examples of pp wave space- 
times are the Aichelburg-Sexl ultraboost [Ij and the Bonnor beam [llj. The 
pp spacetimes are studied not just in the context of general relativity but also 
as purely mathematical constructs in studies of Lorentzian geometry and are 
of interest in string theory due to them often being exactly solvable back- 
grounds. Further they have generated interest due to the causal properties 
of these spacetimes ^J\ [55] . 

As we have seen, purely gravitational plane waves, as opposed to the more 
general pp-wsNes |53| [SH [61], can be characterised by the fact that H{u, x, y) 
is a quadratic function of the transverse coordinates [T2| [22| [61] . 

2.3. Polarisations of the two forms 

We saw in section (I2.1.2p that the general form for a plane gravitational wave 
in the Brinkmann form is 

ds^ = dudv+ {[x^ - y^] H+{u) + 2xyH^{u)] du^ + dx^ + dyl (2.30) 

We see Hj^{u) and Hx{u) are clearly the -|- and x components of the polar- 
isation for the wave. Thus in this form for the metric the two polarisation 
modes are explicitly seen to decouple and to superimpose linearly — quite 
similarly to the situation in Maxwell electromagnetism — and by choosing 
H^{u) and H^{u) appropriately we can construct not just -|- and x polar- 
ized waves, but also circular polarisation, elliptic polarisation, and even more 
general polarisation states. As a simple example, circular polarisation would 
be given by 

ds^ = d-udf -|- H{u) {[x^ — y'^] cos{VLu) + 2xy sm{nu)^ du^ + dx^ -|- dy"^. 

(2.31) 
Equivalently, changing to polar coordinates, in the x-y plane 

ds^ = du dv + r'^ {cos{2(p) H+{u) + sm{2^) H^{u)} dw^ + dr^+r^ d^^. (2.32) 
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So, for instance, to write a circular polarisation we just have 

ds^ = dudv + H{u) r^ {cos(20) cos(r2-u) + sin(20) sin(i7'u) } du^ 

+ dr2 + r2d02. (2.33) 



In contrast the Rosen form cannot be separated in this simple manner, as 
that would require 

Ruu = -\)^ /^ gls - \ /^ 9'bc f"" 9'da] (2.34) 

to easily decouple in some way, and there is no obvious way in which this 
happens. This means there is no clear way to write down arbitrary polari- 
sations in the Rosen form. Even simple, physically important scenarios such 
as circular or elliptical polarisation pose a difficulty. 

This is only a problem for the strong-field case: due to superposition in 
the linearised case we can clearly decouple and add the different polarisations 
to create any arbitrary polarisation. 

Note again that from section (12. 5 p we have the same pattern of non-zero 
components in the Rosen and Brinkmann forms (cf. equations (12. 8p . (12. 9p . 
(I2.10p ). which indicates that some of our intuition regarding polarisation 
modes will also carry over to the Rosen form of the metric. Since these 
two forms are related in this way, and polarisation modes are related to 
the non-zero components of the Riemann tensor, the comparative difficulty 
in constructing polarisations of the Rosen form is somewhat puzzling, and 
suggests that equation (I2.34p should simplify in some way. Investigating this 
curiosity, a simplification of the vacuum equations and hence a way to write 
general polarisations is detailed below. 

2.4. Linear Polarisations 

Consider the relatively simple strong-field gravitational wave metric in the 
-|- linear polarisation [22j. That is, set g^y = so that qab has only two 
nontrivial components, g^x and Qyy, corresponding to oscillations along the x 
and y axes. The resulting metric can be written in the form 

ds^ = dudv + p{u) dx^ + /(m) dyl (2.35) 

The only non-zero component of the Ricci tensor is: 

Ruu = -U-T + —\- (2.36) 
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Although this expression for the Ricci tensor is compact, ultimately this form 
of the metric turns out not to be very tractable. However, if we write the 
metric in the form 

ds^ = dudv + S\u) {e+^(") dx^ + e-^(") dy^} , (2.37) 

then 

Run = -l\^^ ^ + {Xry (2.38) 

This may not initially look very promising, but it is this version of the metric 
that permits us to make the most progress. In particular note that in vacuum 
we have 

X' = 2^~S"/S, (2.39) 

so that 



X{u) = 2 ^/-S"/Sdu. (2.40) 

and the general vacuum wave for + polarisation can be put in the form 



ds' = dudv + S'{u) {exp(2 / y/-S"/Sdu dx 



+ exp ( -2 / ^/-S"/Sdu 1 dr ^ . (2.41) 

Note that as expected from the Brinkmann form we have one free function 
per polarisation mode. 

To take the strong-field gravitational wave metric in the x linear polarisation 
consider 

ds' = dudv + :^-M±^-M[dx2 + dy^] + [f{u) - g\u)]dxdy, (2.42) 

which could also be obtained from equation (I2.35P simply by performing a 
45 degree rotation in the x-y plane. The only non-zero component of the 
Ricci tensor is again: 

Ruu = -ij + -Y (2.43) 

If we now write this metric in the form 
ds^ = dudv + S\u) {cosh(X(M)) [dx^ + dy^] + 2 smh{X {u))dxdy} , (2.44) 
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(equivalent to rotating equation (I2.37P by 45°) then, as for the + mode, we 
have 

Ruu = -I {^^ ^ + {Xry (2.45) 

In vacuum we can again solve for X{u) and now obtain 



ds^ = dudv + S\u) < cosh ( 2 / y/-S"/Sdu J [dx^ + dy^] 

+2sinh ( 2 / ^-S"/Sdu\ dxdyl . (2.46) 

(as above, we could have obtained this form by rotating equation (I2.4ip by 
45°). There is again one freely specifiable function for this x linear polarisa- 
tion mode. 

By rotating the x-y plane through a fixed but arbitrary angle Oq we 
can easily deal with linear polarisation modes along any desired axis. Thus 
there is no difficulty in representing the +, x or any linear polarisation in 
the Rosen form. The difficulty only arises when considering u-dependent 
polarisations. 

2.5. Arbitrary polarisation 

Let us now take an arbitrary, possibly u dependent, polarisation and consider 
the following metric ansatz, which was found through repeated experimen- 
tation: 

ds^ = dudv 

+ S\u) {cosh(S(M))[e+^(«) dx^ + e-^(") dy^] +2smh{B{u)) dxdy] . 

(2.47) 

Note that this reduces to equation (I2.37P when B{u) = 0, corresponding to + 
polarisation, while setting X{u) = (and changing the arbitrary function's 
name from B{u) to X{u)) gives equation (I2.44p corresponding to x polari- 
sation. Furthermore we have sufficient free functions, namely {S{u), B{u), 
X{u)}, to completely saturate the arbitrary 2x2 symmetric matrix gAB{u). 
A brief calculation yields the only nonzero component of the Ricci tensor: 

Run = -U^ y + {B'f + cosh2[i?(n)] {X'f\ . (2.48) 
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Note that B{u) and X{u) have not decoupled — however a hint on how to 
proceed is provided by noting that the 2-nietric, 

dB^ + cosh'^{B)dX\ (2.49) 

is one of many possible ways of representing the metric of the hyperbolic 
plane, H2. 

Let us try a slightly different representation of the same general metric: 

ds^ = dudv + S'^{u) Ucosh{X{u)) + cos{e{u))smh{X{u))]dx^ 

+ 2 sin(^(n)) sinh.{X {u))dxdy 

+ [cosh(X(M)) - cos{e{u)) sm\i{X{u))]dy^} . (2.50) 



Note that setting 6{u) = corresponds to + polarisation, whilst setting 
6{u) = 7r/2 corresponds to x polarisation, and in general 6{u) = Oo/2 cor- 
responds to linear polarisation along axes rotated by an angle Bq. 

Thus this form of the metric introduces a simple physical interpretation 
for the functions, unlike the previous metric. Here, 6{u) describes the angular 
dependence while S{u) can be regarded as an "envelope" function, describing 
the strength of the wave, and X{u) is the shape of the "distortion". Fur- 
thermore we again have sufficient free functions, now {S{u), X{u), 6{u)}, to 
completely saturate the arbitrary 2x2 symmetric matrix gAsiu). 

Now 

Run = ~IU^ + {X'y + smh\X{u)) {e'A , (2.51) 



2 { S 
and the vacuum equations are simply 

4 ^ + (Xf + smh\X{u)) {9'Y = 0. (2.52) 

Let us introduce a dummy function, L[u), and split this into the two equa- 
tions 

4|^ + (L')^ = 0, (2.53) 

and 

(L')' = (X')' + sinh2(X(M)) {e'f. (2.54) 

The first of these, equation (I2.53p . is simply the equation you would have 
to solve for a pure + or x or in fact any linear polarisation. The second of 
these, equation (I2.54p . can be rewritten as 

dL^ = dX2 + sinh2(X)d^^ (2.55) 
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and is just the statement that L can be interpreted as distance in the two 
dimensional hyperbohc plane, H2. 

2.5.1. Algorithm 

This now provides us with a very straightforward algorithm for constructing 
any arbitrary polarisation strong-field gravitational wave in the Rosen form: 

• Pick an arbitrary L{u) and solve: 

4 ^ + {L'f = 0. (2.56) 

• Pick an arbitrary curve in the (X, 6) plane such that L{u) is hyperbolic 
arc- length along that curve: 

dL^ = dX^ + sinh2(X)d^2^ (2.57) 

(alternatively one could consider these two steps in reverse.) 

• This construction then solves the vacuum Einstein equations for the 
metric 



ds^ = dudv + S\u) I [cosh(X(M)) + cos(^(m)) smh{X {u))]dx^ 
+ 2 sin(6'(M)) smh{X {u))dxdy 
+ [cosh(X(u)) - cos{e{u)) smh{X{u))]dy^} . (2.58) 



In this sense we have completely solved arbitrary polarisation strong-field 
gravitational waves in the Rosen form. 

2.5.2. Comparison to electromagnetism 

Note the similarities (and differences) with regard to the familiar Maxwell 
electromagnetism (and with respect to the Brinkmann form). In electromag- 
netism the two independent linear polarisations can be specified by 

E{u) = E,{u)x + Ey{u)y, (2.59) 

with no additional constraints (compare with equation (I2.30p ). Thus an 
electromagnetic wavepacket of arbitrary polarisation can be viewed as an 
arbitrary "walk" in the (E^jEy) plane. 
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We could also go to a magnitude-phase representation {E, 6) where 

E{u) = E{u) cos d{u) X + E{u) sin e{u) y. (2.60) 

(Compare this with equation (I2.32p ). So an electromagnetic wavepacket of 
an arbitrary polarisation can also be viewed as an arbitrary "walk" in the 
{E, 9) plane, where the (-E, 6) plane is provided with the natural Euclidean 
metric 

dL^ = dE^ + E'^de'^. (2.61) 

In contrast, for gravitational waves in the Rosen form we are now dealing 
with an arbitrary "walk" in the hyperbolic plane, H21 rather than in the 
Euclidean plane. Furthermore, because of the nonlinearity of strong-field 
general relativity there is still one remaining differential equation to solve for 
the "envelope", S{u). 

2.6. Circular polarisation 

We can now adopt the above discussion to formulate strong-field circular 
polarisation in the Rosen form. We emphasise again that there is no difficulty 
whatsoever with weak-field, linearised circular polarisation, it is only for 
strong fields that it is difficult to formulate circular polarisation in the Rosen 
form. Circular polarisation corresponds to 

e{u) = fio u] X{u) = Xq. (2.62) 

That is, a fixed distortion Xq with the plane of polarisation advancing linearly 
with retarded time u (with a coordinate choice to set the constant term to 
zero). The envelope function, S{u), must be determined through the Einstein 
equations. Now the metric becomes 



ds^ = dudv + S'^{u) < [cosh(Xo) + cos(r2o u) sinh(Xo)]dx^ 
+ 2 sin(r2ou) sinh(Xo)dxdy 

+ [cosh(Xo) - cos(fio«) smh{Xo)]dy^\ . (2.63) 

The only nontrivial component of the Ricci tensor is then 

Run = -^^ T + sinh2(Xo) Ql\ , (2.64) 



29 



and the vacuum equations are simply 

g-- = - -""''<f)"° g. (2.65) 

To explicitly link this to the algorithmic construction given in section (12.5. ip . 
note that r2osinh(Xo) is the hyperbolic arclength of a circle segment. Now 
we can solve for S{u), 

S{u) = So cos <^ > . (2.66) 

So explicitly, 

d.^ = dudv + S^o (cos 1^^"^™ ^° ^" ' "°^ ^ ^ ' 

< [cosh(Xo) + cos(fio u) sinh(Xo)]dx^ 
+ 2 sin(r2ou) sinh(Xo)dxd?/ 
+ [cosh(Xo) - cos(fioM) sinh(Xo)]dy2 \ . (2.67) 



This now describes a spacetime that has good reason to be called a strong- 
field circularly polarized gravitational wave. 

For consistency, note that the weak-field limit corresponds to Xq ^ 1 so 
that for an arbitrarily long interval in retarded time u we have S ^ Sq, and 
without loss of generality we can set S* ~ 1 . Then 

ds^ ~ dudt> + da;^ + dy'^ 

+Xo I cos(fio u) [dx^ - dy^] + 2 sin(fio u)dxdy\ , (2.68) 

and compare the above (equation (I2.68P ) to equation (ll.57p . 

Further generalisations to elliptic polarisation are tedious, but given the 
significantly more general algorithm of the preceding subsection, reasonably 
straightforward. 

2.7. Decoupling the general form in arbitrary dimensions 

Considering the calculations considered so far, note we have made progress 
with these forms by separating out a unit determinant piece and some enve- 
lope function. Based on this, one might suspect that there is some general 
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decoupling between the overall "envelope" of the gravitational wave and the 
"directions of oscillation". Let us return to considering the metric in general 
Rosen form 

ds^ = dudv + gABiu)dx'^dx^ , (2.69) 

where for generality the x^ , x^ are now any arbitrary number of dimensions 
{d± > 2) transverse to the {u, v) plane. It is easy to check that the only 
non-zero component of the Ricci tensor is still 

Ruu = ~[l /^ 91b - I /^ 9'bc 9^"" 9'nA] ■ (2-70) 

Let us now decompose the d± x d± matrix qab into an "envelope" S{u) and 
a unit determinant metric related to the "direction of oscillation". That is, 
let us take 

9ab{u) = S^{u) gAB{u), (2.71) 

where det(^) = 1. A related discussion can be found in section 109 of Landau 
and Lifshitz [33]. Our goal is to see if we can make the overall "envelope" 
S{u) decouple from gAsiu). To start, note that 

glj, = 2S S' QAB + S^ 9'ab. (2-72) 

and 

g% = 25 S" QAB + 2S' S' qab + 45 S' g^j, + S^ g%- (2.73) 

Therefore 

O O O \ . .O , AB -/ 1 I r^AB s." 



2 ( y + ^^ 1 rf± + 4- [r" 9ab] + [r" g'U (2.74) 



and similarly 



/^ 9'bc ^^^ 9'nA 



g^''[2S S' gBc + S^ 9'bc]9''''['^S S' goA + S^ 9'da 
1 
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^ {4(5 S'fd^ + 45^5'[^^^ g'A^\ + 8%"^^ g'^^ f" 9'bA} 



4('f)'c^^ + 4(|~) \9^''9'AB\^\9^''9'Bcf''9'nA\ (2.75) 



31 



Now combine the results of equations fl2.75p and ( I2.74p 

R — } ^ n^B II 1 ^AB I CD I 

^uu- - <^g 9ab - ^9 9bc 9 9da 

oil QIQI\ Ql 1 

^ + ^)d^-2j[9^'9'As]-\r'9lB] 
+ (I) ' d^ + (I) [9^' 9'ab] + \r' 9'bc f"" 9'ua\ 
= -^-^d^-]p^''9"AB\ 

+ ^[^^^ 9bc f" 9 da] - (f ) [^^^ 9'ab\- (2.76) 

Because we have defined that det(^) = 1, we have as a matrix identity 

{9^'' 9'ab\ = 0, (2.77) 

and differentiating this one more time 

[^^^ 91b] - [^^^ 9bc f" 9'da] = 0. (2.78) 

Therefore 

Run = -d^^-\ [9^'' 9'bc f" 9'ba\ (2.79) 

or more abstractly, 

Run = -d±Y- 1'{^9]-' [g]' [g]-' m . (2.80) 

Note that we have now succeeded in decoupling the determinant (det{g) = 
S'^'^^] effectively the "envelope" S{u) of the gravitational wave) from the unit- 
determinant matrix g{u). This observation is compatible with all the specific 
examples considered above. 

Now, analogously to section (12.5. ip . consider the set SS{lR,d±) of all unit 
determinant real symmetric matrices, and on that set (not a group) consider 
the Riemannian metric, 

dL' = tr{[g]-'d[g][g]-'m}- (2.81) 
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Then 

S"(u) fdlV 




'^•' = --'<^'-^^[tu)>- (2-«2' 



The vacuum Einstein equations then reduce to 



^ = 2y-rfxy; L{u) = j 2^-d^^du. (2.83) 

That is, an arbitrary polarisation vacuum Rosen wave is an arbitrary walk in 
SS{1R, d±), with distance along the walk L{u) being related to the envelope 
function S{u) in the differential relation above. 



To probe the polarisation modes in more detail, a completely analogous but 
slightly more complicated calculation determines the Riemann tensor com- 
ponents to be 



SS'g'^^ . (2.84) 



D I ec"" , o2 ^ '■" ^ -I -CD -/ 

^uAuB — - < o^ gAB + J - 9ab ~ T 9ac 9 9db 

Note that this is compatible with equation 02.51) . 

2.8. Discussion 

The pp waves, and gravitational plane waves in particular, are interesting for 
a number of reasons, including such diverse topics as the study of colliding 
plane waves and quantum theories of gravity. 

In spite of gravitational waves in both the Rosen and the Brinkmann 
forms having been known for many decades they are still of considerable 
interest in research and, as we have found, there are still a number of surprises 
in the Rosen form. 

The two forms that are commonly used for gravitational waves have their 
respective advantages and disadvantages. However, the fact that it is entirely 
trivial to represent arbitrary polarisations in Brinkmann form, but unclear 
how to do so in the Rosen form, hints that the previous understanding of 
polarisations of the Rosen form was incomplete. 

We have found that splitting the metric into an "envelope" function and 
a unit determinant matrix succeeds in decoupling the vacuum equations into 
parts depending only on a walk in some appropriately defined polarisation 
space, and a further differential equation relating this walk to the "envelope" 
function. In (3+1) dimensions this polarisation space is the hyperbolic plane. 
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This decoupling, and the corresponding defining of the polarisation space, 
can be extended to an arbitrary number of dimensions. Furthermore we 
have applied this construction to a case of particular interest, the circularly 
polarized wave. 



r CHAPTER 3 

Penrose limits, ultra-local limits and 
Bianchi cosmologies 



In this chapter we will examine several seemingly unrelated aspects of gen- 
eral relativity. These topics will all be relevant for understanding the next 
chapter, in which knowledge of all the subjects presented will be useful. 

Firstly, we will look at approximations around an arbitrary point and an 
arbitrary null geodesic in any spacetime, the latter of which is known as a 
Penrose limit, and produces the pp wave spacetimes examined in the previous 
chapter. 

Next we will move on to various cosmological models, starting with the 
simple FLRW spacetime, then the set of all homogenous cosmologies. We do 
this not because cosmological concerns are directly important to the research 
presented, but because this is the natural context in which the spacetimes 
useful for the research which follows in the next chapter are most frequently 
encountered. 

We will also look at other sorts of limits: ultra-local limits, which have 
mostly been developed in a quantum gravity context, and some simplifica- 
tions possible in cosmological contexts looking at approximations near sin- 
gularities. 

3.1. The flat space limit 

It is clear that the zeroth term in an expansion for a general spacetime will 
be fiat spacetime, so the limit for an arbitrary spacetime around a general 
point will be a particular fiat space limit. We can precisely describe this limit 
as follows: 

Begin with an arbitrary spacetime 

ds"^ = gabix) dx" dx\ (3.1) 

and implement the limit in stages: 

1. Change coordinates 

x^^Xx", (3.2) 

34 



35 



to obtain 

ds'^ = \'^gabi\x)dx''dxK (3.3) 

2. Perform a confornial deformation 

ds' = — = gai,{Xx)dx''dx'. (3.4) 

3. Now take the limit as A — )■ 

df = gabiO) dx" dx'' . (3.5) 

This construction gives us the flat space approximation to the metric around 
the point x = 0. This could then be used as the zeroth term of an expansion 
to approximate the spacetime close to x = 0. 

Such an expansion is already well-known, using a Riemann normal coordi- 
nate system expansion. Riemann normal coordinates are coordinates where 
at a point (chosen to be af = 0) the metric is r]ab and the Christoffel symbols 
are zero. Such a coordinate choice is always possible. Then the expansion 
about X = is 

gab{x) =r]ab+-RacdbX''x'^+-VeRacdbX''x'^x'' + 0{x^). (3.6) 

3 

For instance, see 

3.2. Penrose limits 

The previous chapter dealt with plane waves. Although some of the applica- 
tions and areas of interest for these spacetimes were mentioned, one context 
in which these spacetimes have been studied was omitted: that they are the 
result of Penrose limits. This limit, suggested by Roger Penrose in 1976 
[54], is the zeroth approximation for a spacetime around an arbitrary null 
geodesic. 

3.2.1. Construction 

To construct the Penrose limit pick an arbitrary null geodesic, 7, in any 
arbitrary spacetime, and in some tube surrounding the geodesic choose a set 
of adapted coordinates, also sometimes referred to as Penrose coordinates, 

d4 = dUdV + a{U, V, f) dV^ + k{U, V, y') dV dy' + g,,{U, V, y') dy' dyK 

(3.7) 
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It is always possible to put a spacetinie into these coordinates locally [S] 
(this is not required globally and as such does not impose any restriction on 
the type of spacetinie considered). This coordinate system corresponds to 
embedding a null geodesic in a congruence of twist-free null geodesies, given 
by V and y* constant. 

Now implement the Penrose limit in stages: 

1. Make a coordinate transform 

y = A%; y' = Xx'; U = u, (3.8) 

then 

ds^ = A^ du dv + X^ a{u, X\, Xy') dv'^ + A^ hi{u, X^v, Xy') dv dy' 

+ X^g,,{u,X^v,Xy')dy'dyK (3.9) 



2. Perform a conformal deformation 

dsi 
A2 

+ g,,{u,x\,Xyl(^y'^y'- (s.io) 



dstf = —J- = dudv + X^ a{u, X^v, Xy') dv^ + 2Xbi{u, X'^v, Xy^) dv dy' 



3. Take the limit A — )■ 0. The net effect is 

ds^ = dudv + gij{u, 0, 0) dx' dx\ (3.11) 

Thus we see that the result of this limit is a pp wave spacetime. This 
could then be used as a zeroth approximation to the spacetime around the 
geodesic in some expansion. 

Therefore this limit takes any lightlike geodesic in any spacetime and 
essentially "blows up" that geodesic to becomes what turns out to be a pp 
wave spacetime. A physical interpretation is that this is the limit of an 
accelerating observer approaching a lightlike geodesic whilst simultaneously 
adjusting their clocks to run progressively faster (ensuring the spacetime does 
not become degenerate), approaching the limit of measuring u [54]. 

Penrose limits have generated interest for a number of reasons. They 
are particularly important in the context of string theory where pp waves 
often provide exactly solvable backgrounds [33| IM] . Some cases that have 
been particularly well-studied are the limits of AdS„ x S" and Schwarzschild 
or more general singularities |51 El El [ID] • 
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3.2.2. Features of the Penrose limit 

An important question for any limit is: what properties of the original metric 
are preserved in the particular limit? These properties are called hereditary 
properties [SI Ell- 
in the case of Penrose limits a property will be hereditary when it remains 
the same through the various steps in the limit (13.81) - (13.111) . In many cases 
we are most interested in coordinate independent properties, and in this case 
we only need to consider whether such properties survive rescaling (I3.10p and 
the limiting process (13. lip . 

One important result is that if there is a tensor field which is zero for all 
A > 0, then this field will also be zero in the limit f3T]. Thus, for instance, a 
Ricci fiat spacetime will have a Ricci flat Penrose limit. 

Penrose limits do not preserve all structures that may be of interest. For 
instance they erase the presence of horizons [55|I5U]. Note that Penrose limits 
do not preserve the form of the stress-energy. We previously noted that pp 
waves have only one non-zero component of the Ricci and Einstein tensors, 
allowing only vacuum or null dust. This means that the Penrose limit of 
any spacetime also only has at most one component of the stress-energy, Tuu 

[HlESlEaEI]. 

3.3. FLRW COSMOLOGY 

Our fundamental understanding of cosmology is based on the evidence that 
the universe is highly homogeneous and isotropic (and has has been back to 
at least the time of the decoupling of the CMB). This gives us the information 
that (averaging over sufficiently large volumes) the matter in the universe is 
well described by a perfect fluid, and 

Rij = i^Rgij, (3-12) 

implying that the curvature of space is constant, leaving essentially only three 
options of fiat, spherical or hyperbolic space. 

The direct consequence of these symmetries is that we have a very simple 
spacetime: the Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime: 

ds^ = -dt^ + a{t)' {y^^ + r' [de^ + sin2(e) d</>2] | , (3.13) 

where the free parameter, a(t), is a scale factor describing the "size" of the 
universe, and 
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• k = 1 ^ positive spatial curvature, spatial slices have spherical geo- 
metry. 

• k = ^ zero spatial curvature, spatial slices are fiat. 

• k = — 1 ^ negative spatial curvature, spatial slices have hyperbolic 
geometry. 

3.3.1. Einstein equations and conservation law 

The Einstein equations produce these relations, describing the evolution of 
the universe in terms of the matter present: 



^7r p = 3 



inp 



It is sometimes useful to rearrange the Einstein equations as: 



3 


'a' k' 

_a2 + a?_ 


■) 


a' 
a" 


' k a 

r + - + 2- 

a"^ a 



(3.14) 
(3.15) 



Rab — Tab ~ 7:T Qa 



(3.16) 



which lets us rewrite the above relations (equations (13.141) and (I3.15P ) as 

4 



and 



3-[p + 3|>] = --, 



dak , , 

- + 2- + 2- = 47rp-p. 
a a^ a^ 



(3.17) 



(3.18) 



Furthermore, if we apply covariant conservation to this spacetime geometry 
we obtain a conservation law, 



dt 



3.3.2. Raychaudhuri equation 



•^ r 31 ^ r 



dt' 



(3.19) 



We can simplify the (timelike) Raychaudhuri equation (as described in sec- 
tion (II. 8p ) by noticing that we have a natural timelike congruence defined by 
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u = dt, while we also have u"" = (1, 0, 0, 0). For this timelike congruence the 
vorticity, shear and acceleration are zero. Furthermore, 

e = -3-. 3.20 

a 

Therefore we can simplify the Raychaudhuri equation to: 

R^hu''u'' = -3-. (3.21) 

a 

In this simple spacetime the Raychaudhuri equation is reduced to a statement 
about the Ricci tensor: 

Rtt = -3 -. (3.22) 

a 



3.4. BlANCHI COSMOLOGIES 

Although the universe we see today is highly homogenous and isotropic an 
interesting question is whether it was always this way, or if it was possible to 
start in a highly inhomogeneous or anisotropic state and then, through the 
action of the Einstein equations and dissipative physics, produce the highly 
isotropic state we see today. Could a large range of initial conditions lead 
to a universe like ours? Another question that made people look outside the 
FLRW metric is the horizon problem, whereby far distant parts of the sky are 
mysteriously in thermal equilibrium. While some of these problems are now 
better understood in terms of an inflationary era, for a time one approach was 
to investigate these issues by using anisotropic (or inhomogeneous) classical 
general relativity, and some of those results are still of considerable interest. 

Due to the need for mathematical simplicity one of the ways this is in- 
vestigated is to initially retain the assumption of homogeneity while allowing 
anisotropy. This may seem like quite a severe physical restriction, but some 
of the results in homogenous spacetimes have since been extended to inho- 
mogenous ones. 

Homogeneous spacetimes are given by the nine types in the Bianchi clas- 
sification [71 |26]. This classification was known before the development of 
general relativity via a study of Lie algebras on homogenous 3D topologies. 
The Bianchi cosmologies are found by adding on a time dependence to the 
3D metric of the classification, and embedding this in four dimensions, so 
that the Bianchi cosmologies can be described by the metric [5^ IF2j 



ds^ = -dt^ + gkiit) {e\ dx') (e'j- dx^) , (3.23) 
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where the basis vectors e" are functions of the spatial variables and are related 
by the Lie bracket 

[efc,Q] = C"^fc,e„, (3.24) 

with C"^ki known as the structure constants which are, by construction, anti- 
symmetric. The Bianchi types are then classified according to the form of 
the structure constants, which can be decomposed as 

C™H = efc;pnP" + 25™[,afc], (3.25) 

with ekip the Levi-Civita symbol, n a diagonal tensor (rii, 712,^3) and aa = 
(a, 0,0). Alternatively, the closely-related objects 

are sometimes referred to as the structure constants. 

The Bianchi types are characterised by whether a is non-zero and the sign 
of the components of n, as shown in table (13. 4p . See [26] for details. 



Type 


a 


111 


112 


ns 


I 














II 





+ 








III 


+ 





+ 


- 


IV 


+ 








+ 


V 


+ 











VIo 





+ 


- 





VI„ 


a 





+ 


- 


VIIo 





+ 


+ 





VII, 


a 





+ 


+ 


VIII 





+ 


+ 


- 


IX 





+ 


+ 


+ 



Table 3.1: The Bianchi classification 



The Bianchi types are also often classified into class A with a = and 
class B with a 7^ 0, so that class A consists of {/, //, VIq, VIIq, VIII, IX}, 
and class B of {III, IV,V,VIa, Vila}- The familiar FLRW cosmology is a 
special case of several of the nine Bianchi types. 
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Furthermore it is possible to write down the Ricci tensor for generic 
Bianchi types as in [50] , 

R% = ~-^9'' - l^g^'^g^'; (3.27) 

" 2 dt ^ 4 dt ^ dt ^ ' ^ ^ 

R'r = -1^9^' (^V - S^C-mj) ; (3.28) 

1 d f n:^9jkui 



With 



R\ = — Jg^g"' - P%. (3.29) 

^ 2^dt V^ dt ^ ' ^ ^ ' 



Tj / /^kl r~i 1^ r^kl r^ (^ kl(~i _, (~ik (~il _, (~ik r~i I 



(3.30) 



3.4.1. Vacuum Bianchi types 



Many studies are done in the case of the vacuum solutions of the Bianchi 
cosmologies, as under many circumstances of cosmological importance the 
overall behaviour is not affected by the form of matter. Looking at this case, 
it is common to use e" = {l{x), m{x), n{x)), so 

ds^ = — dt^ + (a^(t)/j(x)/j(x) + l?(t)mi{x)mj{x) + (? {t)ni{x)nj{x)) dxMa;-'. 

(3.31) 
Then, simplifying the general results, 

-R\ = - + ^ + - = 0; (3.32) 

aba 

Where (A,/i, i^) are the structure constants (G;, C„„, C^m), and the dot in- 
dicates a derivative with respect to t. 



Changing variables can also be helpful; 

a = ln(a); (3 = ln(5); 7 = ln(c), (3.36) 
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with a change of time coordinate 

dr = ^ dt. (3.37) 

abc 

So the system of equations fl3.32p -f l3.35p can now be described by 

2a={i2b^-uc^Y -X^a*; (3.38) 

213 = {Xa^ - uc^Y - 12%*; (3.39) 

27=(Aa2-/i62)'-z/V; (3.40) 

and 

d/3 + d7 + /37 = ^ (Aa^ + fi%* + z/^c^ - 2Xfia%^ - 2Xva^c^ - 2^v}?c^) , 

(3.41) 
where now the dot indicates a derivative with respect to r. 

3.5. BiANCHi Type I 

The Bianchi type I is described by all the structure constants being zero. This 
means that we can always pick a triad ^i such that the three dimensional 
slice is given by Euclidean space. However, we want independence of such 
a coordinate choice when embedding in four dimensions, and so the general 
Bianchi type I is: 

ds^ = -dt^ + g^^{t) dxMx^'. (3.42) 

Note the Bianchi type I includes FLRW with A; = as a special case. 

This is the simplest of the nine forms, and one of the most frequently 
studied (for a selection of classic and recent sources see [131 [HI [HJ [321 [3S1 
[181E21[551[571E9]). 

In the cosnio logical context, it is very common to see the full Bianchi type 
I simplified to a diagonal metric. This is appropriate when one is considering 
the most likely forms of the stress energy of relevance to the large-scale 
universe: dust and perfect fluids (see section I^T|1 . but is less relevant for our 
purposes as discussed below. 

3.b.l. Kasner solution 

The unique vacuum Bianchi type I solution is known as the Kasner metric 

ffl, 

ds2 = _dt2 + ^t2p,^^i^ ^3 43^ 
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with constraints 



J2p^ = ^; Ep' = 1- (3.44) 



Choosing pi < P2 < Pa, the exponents have the ranges 

-l<Pi<0; 0<p2<| 1<P3<1- (3.45) 

Furthermore they can be expressed in terms of the single parameter u: 

-u l + u u + u^ 

This obviously includes the flat three space {pi,P2,P3) = (0,0,1). Other 
than this solution, one of the exponents must always be negative; as t — )■ 
the universe shrinks along two of the primary axes whilst growing along the 
other. Generally this must also happen at different rates: the only other case 
where two exponents are the same is {pi,P2,P3) = (~|, |, §)• 

The Kasner metric has been studied as one of the simplest possible cos- 
mologies, and also in understanding aspects of the approach to singularity in 
more general cosmologies. 

3.6. Ultra-local limits 

Simply put, the ultra- local limit is the first step in a possible strong cou- 
pling or high-temperature expansion for some field theory. The standard, 
most well-developed expansion for field theories is a weak-coupling expansion. 
There is, in addition, a strong-coupling expansion, which has the ultra-local 
limit as its zeroth term. 

This strong coupling expansion has been well examined in the case of 
quantum field theory (see [1[ |30] and chapter 10 of jl2]) and has also been 
investigated as a starting point for attempts at quantising gravity 



Essentially this strong-coupling expansion is done by perturbatively adding 
small gradient terms to an independent evolution of separate spatial points. 
It is thus seen that this can be viewed, from a statistical mechanics perspec- 
tive, as a high-temperature limit. 

To make this mathematically precise, in the case of gravity, take the 
Langragian for the ADM decomposition. 
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The ultra-local limit corresponds to setting the potential term, ^^^R, to 
zero. As the extrinsic curvatures have only derivatives with respect to time, 
all spatial derivatives have dropped out and thus each point has decoupled 
from every other point, and evolves separately. 

We could also explicitly reintroduce units, 

C^N I^V. (}^ K, A-« - 15^ A- + -£^ <»)fl) , (3^48) 

so we see that one way of implementing this ultra-local limit is by taking 
the coupling constant, Gn — )■ oo, making it a strong coupling limit. Alter- 
natively, we could see this as taking the speed of light to zero. This forces 
lightcones to become vertical, so that no two points can interact. 

3.7. Approach to singularity 

Bianchi cosmologies, and type I cosmologies in particular, have often been 
of interest in examining the approach to a general singularity. The general 
idea is that close to the singularity the zeroth approximation to a general 
spacetime is a Kasner metric, with spatial curvature terms added on in some 
simple way. There are two ways this could happen; a smooth approach to 
the singularity (velocity-dominated) or chaotic (Mixmaster behaviour). In 
general, close to classical singularities, it is believed that either 



• The cosmology is Kasner-like. 

• The cosmology is oscillatory, "generalized Mixmaster" behaviour. 

• Something else. 

With many believing that generically the behaviour is Mixmaster. There has 
been a significant amount of work into discovering when which behaviour 
happens and if and when the "something else" might occur. 

3.7.1. Velocity dominated limits 

A velocity dominated cosmology, (which is also referred to as having a ve- 
locity dominated singularity, or being an asymptotically velocity-term dom- 
inated (AVTD) cosmology), is a cosmology that can be approximated on a 
coordinate patch by a simpler, ultra-local cosmology near the initial singu- 
larity j231 SH SS] • Denote an ultra-local metric by cjab and extrinsic curvature 
by Kab- A cosmology, Qab with extrinsic curvature Kab is velocity dominated 
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if, when put into synchronous coordinates, in some precise way (for details 
see [231 E]), as t ^ 0, 

9ab -^ 9ab; Kab -^ Kab, (3.49) 

This imphes that close to the singularity the potential terms can be ig- 
nored, so the evolution can be approximated by the kinetic terms |28j. 

Note that this approximation changes the form of matter under consid- 
eration, discarding terms in the stress-energy to permit one to approximate 
by the ultra-local metric. 

Of particular interest is the case when the cosmology is approximated by a 
Bianchi type I metric near the singularity. These cosmologies are sometimes 
referred to as Kasner-like, or as exhibiting Kasner-like behaviour [21j. A 
large number of spacetimes are, or are conjectured to be, Kasner-like [6j. 

A possible source of confusion is that sometimes Mixmaster behaviour 
(see below) is referred to as being velocity dominated [SB], though many 
researchers regard these as competing behaviours |2H l66]. 

3.7.2. BKL singularities and Mixmaster dynamics 

A cosmological concept related to these ideas has emerged from examinations 
of the vacuum Bianchi type IX spacetime (type VIII is often also included in 
these calculations as the behaviour is very similar), examined independently 
by Misner [19] and Belinsky, Khalatnikov and Lifshitz [3]. 

For vacuum Mixmaster, we use X = i' = fi = lm equations (I3.27p - (l3.29p . 
That is: 



2a={b'-c'y-a'' 
2^={a'-b'y-c' 



(3.50) 
(3.51) 
(3.52) 



a$ + aj + $^ = j (a^ + fe^ ^ c^ _ 2a%'^ - 2a^c^ - 2\?^) . (3.53) 

If there is some point in time when the righthand side of equations O3.50p - 



(I3.53P could be neglected, Kasner behaviour would occur, so that 

aoct^^; fcoct^^ c(xt^\ (3.54) 

However, the righthand side will act as a perturbation, so the spacetime 
would not stay in this Kasner regime forever. We know from an examination 
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of the Kasner spacetinies (section f l3.5.ip ) that in this case as t — > two of 
the terms would decrease (contracting to the big bang) and the other would 
increase, so the effect of the perturbation will act through the term that 
increases. Arbitrarily pick this to be a. Now, 

a = -\e'^- /3 = 7 = ^e^", (3.55) 



which can be explicitly integrated: 

2 _ 2\p,\n 

"" "cosh(2bi|l]r)' ^'^■^^' 

fe2 = 6^exp(2fi(p2- bi|r))cosh(2|pi|fi); (3.57) 

c^ = cgexp(2n (p3 - |pi|r)) cosh(2|pi|n). (3.58) 

If we are only interested in the behaviour close to the singularity, then we 
can approximate these to 

a (X exp(— fipir); (3.59) 

6ocexp(n(p2 + 2pi)r); (3.60) 

c oc exp(n(p3 + 2pi)r); (3.61) 

with 

t ocexp(fi(l + 2pi)r). (3.62) 

So it is possible to say that 

aoctP'^; feoct^™; c oc t^", (3.63) 

for all times (close to the singularity) with an explicit mapping between the 
primed and unprimed exponents. For details see [3l |50]. Thus we see that 
perturbing a Kasner mode leads to another, different Kasner mode, allowing 
the simplification of a Mixmaster universe into a series of these Kasner modes. 
This mapping is such that the exponents act like Kasner exponents for signif- 
icant lengths of time, with sudden changes to new exponents. Furthermore 
there are generically an infinite number of these oscillations as t — )■ 0. 

It was conjectured j3] that this behaviour describes the approach to singu- 
larity not just for vacuum Bianchi type IX, but for generic homogeneous and 
inhomogenous spacetinies containing most relevant forms of matter. Much 
investigation has gone into trying to determine how true this conjecture is. 
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Essentially it is argued that the Mixmaster dynamics of Bianchi IX space- 
times characterise a generic classical big bang. This is the idea of the BKL 
singularity. 

In the case of inhomogeneous cosmologies, it is believed that generically 
the approach to singularity will demonstrate Mixmaster behaviour in the 
sense that different spatial points will approach different Mixmaster solu- 
tions [66j. However, few investigations have looked at whether Mixmaster 
behaviour occurs in inhomogenous cosmologies. 

It is important to keep in mind that these arguments were developed 
heuristically, and as such a number of opposing ideas have developed over 
the many years this issue has been investigated. A collection of useful recent 
articles regarding BKL singularities and mixmaster dynamics are [2T1 |36l EOj 

Eg. 

3.7.3. Kinetic approximations 

We see from the above that it has often been desirable to approximate a 
general spacetime by just the kinetic terms in some appropriate way. Un- 
der certain conditions a generic spacetime can be approximated by a simpler 
one, containing only kinetic terms, after which spatial gradient terms could 
be added on perturbatively. In the case of velocity dominated limits, the 
potential terms can be thrown away close to the singularity, while with Mix- 
master behaviour the potential terms are important only over brief intervals 
where they act to change between Kasner modes. 

3.8. Discussion 

We see that several limits have been developed which involve some form of 
ultra-locality, in which Bianchi type I spacetimes have frequently been of 
central importance. Furthermore taking a limit along a null geodesic has 
generated significant interest. In the next chapter we combine these two 
developments to construct an ultra-local limit along a timelike geodesic. 



r CHAPTER 4 

Any spacetime 

has a Bianchi type I cosmology 

as a limit 



In this chapter we will construct a particular ultra-local limit along an ar- 
bitrary tinielike geodesic, research that was presented in |i20j. The limit we 
describe is closely related to the Penrose limits, and we shall explore the 
properties of this limit and its resulting spacetime. 

4.1. Construction of the limit 

To define the ultra-local limit pick an arbitrary timelike geodesic, 7, and, 
in some timelike tube surrounding the chosen curve, choose a specific set of 
adapted coordinates coordinates (t, x*) such that: 

• The curve in question lies at x* = 0. 

• At X* = the coordinate t is just proper time along the curve. 

• The coordinate system is synchronous. 
This coordinate choice is enough to set 

ds^ = -c^ dt^ + gij{t, x) dxMx^ (4.1) 

where we have explicitly introduced the speed of light. Note that this co- 
ordinate system is only required locally, not globally, and hence this choice 
does not impose any restriction on the type of spacetime being considered. 
Similarly to section (13.2. ip this now corresponds to embedding the geodesic 
into a congruence of twist-free timelike geodesies. 

Now implement the ultra-local limit in stages: 

1. Make the coordinate transformation x* — )■ ex*. Then 

ds^ -^ -c^df + e^g,j{t, ex) dxMx^'. (4.2) 
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2. Change the speed of hght c — )■ ec. So the net effect is 
ds^ ^- -e^c^dt^ + e^5fy(t, ex) dxMx^'. 



(4.3) 



3. Perform a confornial transformation ds^ — t- e ^ds^. So overall 

ds^ -> -cMt^ + ^i^(t, ex) dxMx^'. (4.4) 



4. Take the limit e — )■ 0. Then the net effect is 

-.2 > „2j+2 



ds' -> -cMr + c/ij(t) dxMx^. 



(4.5) 



This definition is clearly "ultra-local" — one has effectively set c — )■ so 
that nothing propagates, such that each point in "space" becomes discon- 
nected from every other point and evolves separately; and each individual 
point in space is then "blown up" by the confornial transformation to yield 
an independent universe. Essentially we have taken this limit of the space- 
time to be given by the geodesic itself, isolated from the influence of the rest 
of the spacetime. 

An alternative (but completely equivalent) picture of our limit can be con- 
structed by again starting from 



ds = — c dt + gij{t,x) dx^dx^ . 



(4.6) 



Consider the decomposition of Qat, as in section (II. 8p . into a projection met- 
ric, hab, and timelike vectors Ua- 



gabit, X) = habit, X) - UaUb, 



with 



h 



ab 



' 





_ 


9ij _ 



M„ = (c, 0,0,0). 



Now consider the three dimensional slice 

dg2 = habit, x)dx''dx'' 
and again implement the limit in stages 



(4.7) 
(4.8) 

(4.9) 
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1. Take the coordinate transformation 

t^t; x' -^ ex\ (4.10) 

resulting in 

dg2 = e'^habit, ef)dx"dx^ (4.11) 

2. Perform a conformal transformation 



dg,2 = ^ = habit, ex)dx''dx\ (4.12) 



3. Now considering 
take the hmit e — )■ 



ds^ = -cMt^ + dg,^ (4.13) 

dsl = -c^df + Kbit, 0)dx'^da;^ (4.14) 



or 



ds^ = -c^dt' + gij{t, 0)dxMx^. (4.15) 



From this point of view we have decomposed the spacetime and then taken 
the flat space limit as in section (13. ip on the three-dimensional slice only. In 
this coordinate system, we are left with the value of the metric along the 
timelike geodesic. 

Note that the output from this ultra-local limit is exactly the definition of a 
general Bianchi type I spacetime — the general definition before you make 
any assumptions about diagonalisability of the spatial metric gij{t). Also 
note that if one starts with a congruence of timelike geodesies in the original 
spacetime, then the output of this ultra-local limit is a collection of general 
Bianchi type I spacetimes, one Bianchi spacetime being attached to each 
timelike geodesic. 

4.2. Features and relation to other limits 

Clearly this limit closely parallels the Penrose limit construction, in fact one 
can argue that it is as close as one can possibly get to a "timelike Penrose 
limit" — though one should certainly not expect all of the features of the 
usual "lightlike Penrose limit" to survive in this timelike case. 
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Note carefully that this ultra-local limit does not share the hereditary 
properties discussed in section (I3.2.2P of the ordinary Penrose limit. The 
information about the original spacetime retained or discarded in the two 
limits will be different. In particular, note the ultra-local limit of a Ricci flat 
spacetime is not necessarily Ricci flat — from the ADM viewpoint (described 
in section II. 9p this is with hindsight obvious since during the limiting pro- 
cess the extrinsic curvatures and the intrinsic (spatial) curvatures scale in a 
different manner. 

In terms of the physical usefulness of this limit, it is best looked at in the 
context of two other well-known limits: 

• If we are interested in approximating a spacetime at a point the ap- 
propriate construction would be to choose Riemann normal coordinates 
and then approximate the spacetime by the tangent space at that point, 
as described in section fl3.1D. 



• If one is interested in approximating a spacetime along a null geodesic 
then the appropriate construction is to use Penrose's adapted coordi- 
nates and take the Penrose limit along that null geodesic as described 
in section (13. 2p . 

• Likewise, if one wishes to approximate a spacetime along a timelike 
geodesic, the appropriate choice is to use synchronous coordinates and 
then take the limit described above along the timelike geodesic as de- 
scribed in the process above. 

4.3. General Bianchi type I spacetime 

Now that the speed of light c has done its job, let us again adopt units where 
c = 1, so that we are considering 

ds^ = -dt^ + gij{t) dxMx^'. (4.16) 



Geometrically this is a rather simple spacetime. Define the extrinsic curva- 
ture 

K^J = -i%, (4.17) 



then in particular 

tr(K) = - 

V93 dt 



tr(i^) = -^ ^, (4.18) 
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where the notion of "trace" imphcitly involves appropriate factors of the 3- 
nietric and its inverse. We can speciahse the usual results for the ADM 
decomposition (for an explicit computation, see (and appropriately modify) 
the discussion of the ADM decomposition in Misner, Thorne, and Wheeler 
|47] ) . to determine the useful geometric tensors for the Bianchi type I space- 
time. 

4.3.1. Riemann tensor 

The key results for the Riemann tensor are 

Rijki = KikKji - KaKjk, (4.19) 

and 

with all other components zero. 

4.3.2. Ricci tensor and scalar 



A brief computation yields 



while 



Therefore 



i?oo = ^^-tT{K'), (4.21) 



dK 

R^J = -^-2K,^K'^,+tT{K)K,^. (4.22) 



R = -2^^^+[tT{K)]^ + tT{K^). (4.23) 

iJiL 



4.3.3. Einstein tensor 



The Einstein tensor can also be evaluated from the above, or can be read 
off from Misner, Thorne, and Wheeler [4^ (see (21.162 a,b,c) p 552). The 
time-time component is simple 

G,o = I {MK)r - tT{K')} . (4.24) 
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In contrast the space-space components are relatively messy 
Gr, = - ^ [K^j - g^MK)] + 3tr(K)K,, - 2{K^ 



itj 



lg,,{MK)r + tT{K')}. (4.25) 



4.4. Simplifying the general Bianchi type 1 spacetime 

These calculations for the relevant tensors of the general Bianchi type I space- 
time are reasonably simple but there are several other tricks that may be 
employed to put them into forms useful for further exploration of these space- 
times. 



4.4.1. Mixed components 

It is often advantageous to work with mixed components such as R^h. Work- 
ing from the above, (or suitably adapting results from the key article on BKL 
spacetimes [3]), one has: 

R% = -^^ + tTiK'y, (4.26) 

R\ = 0; R^ = 0; (4.27) 



R'.—^}M^K'.]- (4-28) 



Also note from equation (I4.28P and using equation (I4.18P that (summing only 
over the space indices) 

R'. ^ -i=^v^. (4.29, 



For the Einstein tensor, the space-space component (equation (I4.25P ) simpli- 
fies somewhat, becoming 

r^ 1 ^ I /TTK^]^;^^ IM^ ^^(-^)]' t^(^')\ rA-^m 
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4.4.2. Unimodular decomposition 

We will now attempt a decomposition similar to that in section (12. 7p : 

gij{t) = a{tY g,j{t), (4.31) 

where det(^jj) = 1. Here a{t) can be viewed as an overall scale factor (similar 
to that occurring in FLRW cosmologies), while cjij describes the "shape" of 
space. Then 

^ = 2aag,, + a'^, (4.32) 
and so from equation (I4.17p . 

Kij = —ad Qij + c? Kij, (4.33) 

with g"^^ Kij = because of the unit determinant condition. Now define 

K',=g"'kk, (4.34) 



so that 



Thus 



Then 



and 



K'. = --6% + K'.. (4.35) 

a 



tT(K) = -3 -; (4.36) 

a 

tT{K^) = 3 d^a^ + tT{k^). (4.37) 



Roo = 3- + tT{K^); (4.38) 

(Jj 



Rij = a^ I -^ - 2krkmj + tiik)kij \ + {ad + 2d^)gij, (4.39) 

whilst the Einstein tensor is given by 

Gij = a?Gij + aakij - (2ad + d^)gij. (4.40) 

We find these do not simplify as well as the expression in equation (I2.79p . 
as we no longer have a trace over all indices. However it is possible to make 
a slight improvement by combining this decomposition with the previous use 
of mixed components. Then 

i?% = 3- + tr(i^2); (4.41) 

a 



55 

and 





'k^, 


J a [a a? } 


(4.42) 


Summing over the space components only we have (in agreement with the 


previous subsection, see equation (I4.29P ) 




-■-{^4}^^T■ 


(4.43) 


For the Einstein tensor the interesting piece is 






K'j_ 


a -^ 2 ya a^ ] ■' 


(4.44) 


4.4.3. Summary 





We have attempted to find a simple way to write down various useful ten- 
sors for the general Bianchi type I spacetime. The spacetime curvature of 
the general Bianchi type I spacetime can be evaluated as simple algebraic 
combinations of the 3x3 matrices gij , ijij , and cjij . If we split the geometry 
into an overall scale factor a{t) and a shape Qij, then the curvature can be 
evaluated as a simple algebraic combination of a, a, a and the 3x3 matrices 

9ij) 9iji s-iid 9ij- 

Due both to ease of physical interpretation and simplicity of equations it 
is often useful to work with mixed components, and to split the metric into 
a scale factor and a unit determinant piece. 

4.5. Stress-energy in the ultra-local limit 

What happens to the stress-energy in the ultra-local limit we are interested 
in? Based on physical intuition one expects that for the energy-momentum 
flux Toi — )■ 0, as the limit means that each point is disconnected from all 
others, evolving independently, but what about Too and Tjj? Remember that 
for the light-like Penrose limit the stress-energy also changes as Tab -^ T^u, or 
alternatively Tab ~^ ^a^6 with k some null vector, with all other components 
vanishing |S1 ESI [Ml E] • 

It is convenient to develop the ultra-local limit in a somewhat different 
form: 

1. Make the replacement 

gabit, x) -^ gl^{t, x) = Qabit, ex). (4.45) 
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2. Furthermore, for any generic field \l/(t, x) make the replacement 

^(t, x) -^ ^e(i, x) = ^(t, ex). (4.46) 

3. Then consider the limit as e — ?■ 0. 

We consider a few specific cases. Most of the forms of stress-energy de- 
scribed below were considered in section (11. 5p . 

4.5.1. Perfect Suid 

Consider a generic perfect fluid in a generic spacetime, 

T"'' = (p + p)V^V'' + pg^\ (4.47) 

and ask what the ultra-local limit might be? Working from the geometrical 
side we have G"^* — )■ in the naturally adapted coordinate system. Thus via 
the Einstein equations it follows that we must also have T°* — )■ 0, which in 
turn implies \^* — )■ 0. Therefore 

Too^p; Toi^O; T^j ^ p g^,. (4.48) 

That is, the generic perfect fluid reduces in the ultra-local limit to a co- 
moving perfect fluid. 

4.5.2. Scalar field 

Consider a scalar field 0(t,x), which has stress-energy tensor 

Tab = (jy,a<f>,b - Igab { ( V0)' + V{<j)) } . (4.49) 

Now take the ultra- local limit. Note 

= ^{<P' + e^g''d,4>dj^ + v] 

^^{(/.2(t,0) + \/}, (4.50) 



whence ultimately 



T^o^U<l^''it,0) + V\. (4.51) 



Similarly 



Furthermore 
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T^i = \ {kd.ct>?^ ^ \ { ect>d,ct> \ ^ 0. (4.52) 



^ e29,09,0 - \ {(-02 + e^dct^f) + V] , (4.53) 

whence ultimately 

T^.^[^\t,Q)-v]g,,. (4.54) 

The key point is that the structure of the stress-energy tensor guarantees 
that in the ultra-local limit 

Too -^ P; Toi -)■ 0; Tij -> p gij. (4.55) 

Similarly to the limit of the generic perfect fluid, the scalar field situation 
results in a co-moving perfect fluid. Such simple behaviour is common, but 
by no means universal. 

4.5.3. Electromagnetic field 

For the electromagnetic field the most natural form of the ultra-local limit is 
to consider the vector potential of equation (ll.30p and take 



A"(t, x) -> A^(t, x) = U(t, ex); A{t, ex)) , (4.56) 

in which case 

E^~A{t,Oy, B^O. (4.57) 

So for the stress-energy tensor, 

roo^^l^'^Aj; (4-58) 

Toi ^ 0; (4.59) 

T,. ^ E,E, - ^ {g^'EkEi}g,,j. (4.60) 

In terms of the electric field E, and the naturally defined inner product 
g{E, E), it is best to summarise this as: 

Too^^g{E,Ey, (4.61) 
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Toi -^ 0; (4.62) 

T,,-^E,E^-]^g{E,E)g,j. (4.63) 

Note that in general the 3x3 matrix of space-space components of the stress- 
energy will not be diagonal. Thus this gives an important case where many 
of the cosmological calculations cannot be blindly followed. 

Note further that in reference [59j a somewhat different limit is taken, 
because those authors are interested in different physics. This is done because 
if one is interpreting the Bianchi type I spacetime as a cosmology, as an 
approximation to the large-scale structure of our own physical universe, then 
for physical reasons one might wish to allow i? 7^ 0, while expecting that 
£" — )■ 0, as magnetic fields may be important in cosmology but electric fields 
almost certainly are not. Nevertheless, even in this very different context, 
one key point remains: The 3x3 matrix of space-space components of the 
stress-energy need not and generally will not be diagonal. 

4.5.4. Anisotropic fluid 
Consider a generic anisotropic fluid 

T"^ = (p + pjVV^ + pg""^ + e''^ (4.64) 

Here B is the anisotropic stress tensor, which is always 4-traceless and 4- 
orthogonal to the 4- velocity V . What might the ultra-local limit be? However 
we choose to define this limit, since we know from the geometrical side that 
G"^* — 7- 0, we must via the Einstein equations have T°* — > 0, so that the 
stress-energy tensor is (3+1) block diagonal. This in turn implies that we 
can without loss of generality redefine our variables so that l^* — )■ 0, while 
V^ — )■ 1, and simultaneously choose 0°" — )■ along with 0°* — )■ 0. Therefore 

Too -^ P; T^i -> 0; Tij -> ttjj = pgij + 9^^-. (4.65) 

Once again note we do not necessarily have a diagonal metric for more general 
forms of matter. 

4.5.5. Stress-energy conservation 
Consider the covariant conservation law 

V.T"" = ^ -^da{V^4 T"') + r^.T^'^ = 0. (4.66) 
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For Bianchi type I the only nontrivial component is 



which imphes 



that is 



^ a,(V^r°0) + rO,,T^'^ = 0, (4.67) 



^ 5t(v^p) + ^%vr^^' = 0, (4.68) 



^ ^^^^^"^ 2 dt 



^^.(v/^p) - -l^^^^. (4.69) 

To see the connection with the well understood FLRW spacetime (section 
13. 3p note Qij = a^ gij and then 



- dt{a'p) = -- [2aag,, " «^ ) ^''' (^-^O) 



which we can write as 



i2a„ ^ij ^ ^5'^9ij ij 
1 n n ■ ■rr -J — — n •tt •' 



dt{a'p) = -a'dg.^n^^ - ^a'^n^^. (4.71) 

That is, (defining ttj^ = a^TTij, so that 7r*j = 7r*j and vr*-' = a^'^if'^^), 

^ = -3a>-la3^{*.-p<)"}. (4.72) 

Finally, (defining 6jj = a^ Qij, so that 6*^ = &j and B'-' = a~^ 6*-^), 

dt ^ dt 2 dt ^ ^ 

The first 2 terms are the familiar FLRW conservation law as given in equation 

(I3.19p . The last term involves the trace-free anisotropic part of the (spatial) 

stress tensor, together with the trace-free contribution from dcjij/dt. We can 

also write this as 

dp ^ , , da 

So the anisotropy modifies the conservation equation, splitting essentially 
into parts concerning the change in density due to the change in "size" and 
that due to a traceless contribution due only to the change in shape of the 
space. 



-3{p + p)— + k,^Q'^. (4.74) 
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4.5.6. Summary 

In general the message is that in the ultra-local limit the stress-energy satisfies 

Too -^ P; Toi -> 0; Tij -> ttj^; (4.75) 

and that one cannot say anything more than this without specifying the 
particular form of the matter content. It seems that looking at stress-energy, 
we obtain a more complex picture than the lightlike Penrose limit, whilst 
sharing the fact that both limits erase some terms from the stress-energy. 
Furthermore in this ultra-local limit there is a natural generalisation of the 
energy conservation law normally applied to FLRW spacetimes, with an extra 
term coming from the interplay between anisotropics in the stress tensor and 
changes in the shape (not volume) of the spatial slices. 

4.6. Einstein equations 

In the ultra-local limit the Einstein equations reduce to three significant 
pieces of information — arising from the the time-time component, the spatial 
trace, and the anisotropic part of the spatial trace. As in section (13.3. ip . it 
is most convenient to rewrite the Einstein equations as 



in which case 



Rab 


= SirGNiTab - 


- 2^.- 


-} 




Roo = 4:7r{p + 


3p); 






R\ = 12n{p - 


-P); 




Rij "3-^ kOij = 


STrQij^ 




ion 


(|477)) implies 







(4.76) 

(4.77) 
(4.78) 

(4.79) 



- = -^{p + 3p)-\tTiK'). (4.80) 

a 6 3 

That is, as compared to FLRW (cf equation I3.17P spacetimes, changes in the 
shape of the spatial slices can now contribute to deceleration. The second 
equation (I4.78P yields 

- + 2^ = 47r(p-p). (4.81) 
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So this combination of terms is what one might expect anyway for any spa- 
tially flat FLRW cosmology (cf. equation 13. ISp . Finally we obtain an addi- 
tional relation from equation ( I4.79p . 



dt 



K'j -3-K'j = 87re'j. (4.82) 



Alternatively, 

--^ a^K'j =87ra^e'j. (4.83) 

That is, evolution of the shape of the spatial slices is driven by anisotropics in 
the stress tensor. This is the simplest and most straightforward decoupling 
of the Einstein equations we have managed to find. 

4.7. Conditions for a [block]-diagonal metric 

A key issue in taking this ultra-local limit (and in Bianchi type I cosmologies 
generally) is the question of when it is possible (or desirable) to restrict atten- 
tion to diagonal spatial metrics (or more generally block diagonal metrics). 
In the cosmological setting it is common to see the full Bianchi type I uni- 
verse quickly reduced to the diagonal case. Fortunately, within the context 
of Bianchi type I spacetimes it is possible to prove the following qualitative 
results: 

(diagonal stress-energy) -^ (diagonal metric). 

(block diagonal stress-energy) -^ (block diagonal metric). 

Proof (<^): A [block]-diagonal metric gij{t) implies that gij{t) is [block]- 
diagonal. This in turn implies that Kij{t) and Kij{t) are [blockj-diagonal. 
Through the calculations for the Ricci tensor, this implies Rij{t) is [blockj- 
diagonal, and this then implies Gij{t) is [blockj-diagonal, which implies TTij{t) 
is [block] -diagonal, implies 7r*j(t) is [blockj-diagonal. (Completely equiva- 
lently, this implies Ojj(t) is [blockj-diagonal, which implies 6*j(t) is [blockj- 
diagonal) . 

Proof (^): At time to there is no loss of generality in picking coordinates 
such that both [go]ij = Sij, and such that [go]ij is diagonal. Then without 
loss of generality there is a coordinate system such that [Kojij and [-K'ol^j are 
diagonal. 
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Our definition of [block]-diagonal stress-energy will be to assert that in 
this particular coordinate system T*j(t) is [block] -diagonal, whence G*j(t) is 
[block]-diagonal, and so R^j{t) is [block]-diagonal. 

Note that we have only imposed a restriction on the metric at one par- 
ticular time. We show in this case the Einstein equations do not lead to 
evolution away from a diagonal metric. 

Now, using 



which we derived earlier (equation (14.281) ). as K^j is [block]-diagonal we re- 
arrange to 

^ ' ^K'j] oc {[block]-diagonal(t)}*j-, (4.85) 



dt' 
and so 

^sK'j = ,/9^[KoY, + {[block]-diagonal(t)}V (4.86) 

By our initial assumption on the diagonal nature of [-ft"o]*j we have 

K'j = {[block]-diagonal(t)}^j. (4.87) 

That is 

Vt : g'"" gmj = {[block]-diagonal(t)}V (4.88) 

but then 

Vt : gij = gim {[block]-diagonal(t)}'"„ (4.89) 

which can be formally integrated in terms of a time-ordered product 

[9it)]ij = [f7o]fmexp[{[block]-diagonal(t)}]"^, (4.90) 

= (5imexp [{[blockj-diagona^t)}]™" j, (4.91) 

where in equation (I4.9ip we are again using coordinate choice that the initial 
metric is Sim- Now 

[git)]i, = exp[{[block]-diagonal(t)}].^. (4.92) 

= {[block]-diagonal(t)}ij, (4.93) 

implying that [g{t)]ij remains [block]-diagonal for all time. 
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To summarise: a necessary and sufficient condition for a [block]-diagonal 
metric (in the context of Bianchi type I spacetimes) is that, in the coordinate 
system where [go]ij = Sij and [Ko\ij is diagonal, we have 

7r*j(t) is [block]-diagonal. (4.94) 

The metric and stress tensor then decompose into direct sums: 

9ij = ©A [hAlif, T^ij = ®A {pA[hA]ij}; Tr'j = ®A {pa[SaYj}- (4.95) 
4.7.1. Examples 

For a fully diagonal metric in three space dimensions the three canonical 
examples of this behaviour are: 



1. Perfect fluid: 
That is: 



TTi 



P9ij- 



TT 



p 


" 


p 








p _ 



pS'r 



(4.96) 



(4.97) 



The special sub-case p = corresponds to dust, and the even more 
special sub-case p = p = corresponds to vacuum (Kasner solutions) 
— with all of these special cases (perfect fluid, dust, vacuum) being 
very well studied in a cosmological setting. 

2. Electromagnetic field with the electric field E being an eigenvector of 
9- 



TT, 



V 



Cl9^J + C2V,Vf, V'^V,. 



This is equivalent to: 



TTi 



That is: 



Pi 9ii 

P2 922 

P2 933 



n 



9ij 



pi 

P2 
P2 



911 

922 
933 



(4.98) 



(4.99) 



(4.100) 
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3. General diagonal metric: 

■Kij = diag{pi,p2,P3); gij = diag{gu,g22,g33)- 
This is equivalent to: 



TTi 



Pig 11 

P2 fi'22 

P3 5-33 



That is: 



9ij 



TT" 



Pi 

P2 
P3 



(4.101) 



^11 
g22 

(733 



(4.102) 



(4.103) 



There is also a partially diagonalisable case corresponding to one irre- 
ducible 2x2 block plus a singleton 1x1 block, and the fully generic non- 
diagonalizable case corresponding to the stress tensor being a single irre- 
ducible 3x3 block. These are the only relevant cases for 3-diniensions. 

4.8. Raychaudhuri equation 

In a Bianchi type I spacetinie we can simplify the standard (timelike) Ray- 
chaudhuri equation [Mj [65j by noticing that we have a natural timelike con- 
gruence defined hy u = dt, while we also have u = — (dt)", or more prosaically 
m" = (1, 0, 0, 0). For this timelike congruence the vorticity is zero, w = 0, as 
is the acceleration ^ = u^VhU"" = 0. 

Therefore in terms of the expansion tensor 6ab = ^{aUb), (not to be con- 
fused with the anisotropic stress Qab), we have 



ab 



RabU'^u" = -OabO^' + 



d^ 

dt' 



(4.104) 



Separating out the unit determinant metric, gij = a'^cjij, and using the defi- 
nitions of shear and expansion given in section 01. 8p we find 



Oab = ad gab + a^Kab] 
o'ab = a Kab] 



(4.105) 
(4.106) 
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a^ = ltT{k^). (4.107) 



2 



Furthermore K = , and hence 

a 



, d d^ a a^ 

a dr a a^ 



= 3-; — = 3--3— . 4.108 



So the Raychaudhuri equation is given by 

RabU^u^ = ~tT{K^) - 3 -. (4.109) 

Compare this to both the result for the FLRW cosmology, previously given 
in equation (I3.2ip . 



a 



RabU^u" = -?>-, (4.110) 



a 



to which the Bianchi type I spacetime reduces when Kab = 0, and to the 
result stated by Collins and Ellis |16] for perfect fluid Bianchi cosmologies, 



^(p + 3p) = -3^-2a2. (4.111) 

Attempting to put the Raychaudhuri scalar into a "nice" form, the best pos- 
sible form seems to be 

Rabu'^v!' = i (p + 3p) + ii{k^), (4.112) 

or alternatively 

Rabu''v' = -Ap + ?>p)-\-. (4.113) 

4t Zi Qj 



4.9. Some examples 

We will now look at a few simple examples of this limit. 
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4.9.1. FLRW 

Let us look at the simplest example of this limit that is not completely trivial. 
Consider the FLRW spacetime in Cartesian coordinates, 

This is already in the desired adapted, synchronous coordinate system. The 
only effect of the ultra local limit is that 

1 + i<^:lif±i!) ^ 1, (4.115) 

and so the limit is simply the spatially flat {k = 0) FLRW spacetime. 

4.9.2. Circular orbits in spherical symmetry 

If we have spherical symmetry, the coordinate system adapted to a circling 
observer will have a time independent gij, as throughout the path of the 
geodesic there is no variation in the gravitational field. 

Now taking the limit in the appropriate manner, we remove the depen- 
dence on the spatial variables. Since our final Bianchi type I spacetime has 
no time dependence, we can simply rescale coordinates to get the standard 
form for the Minkowski spacetime. 

4.9.3. Radial infall for the Schwarzschild black hole 

To take the ultra-local limit in the Schwarzschild solution, we first need the 
metric in synchronous form. Let us start with the standard/Schwarzschild 
form: 

ds^ = - (l - ^) dt^ + ^wd^' + ^' (d^' + sin'(^) d^^) . (4.116) 

Now, closely following both section 102 in [43j, and [46j, perform a coordinate 
transform of the form 



t-p= 7 ,.,, , dr, (4.117) 

which defines p. Furthermore define a new time coordinate, r, by 

T-p= I ^dr. (4.118) 

J (2Mr)^ 
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So that 



dr = (dr — dp ) 



r{r,p)={-{T-p)] (2M)3. 



We see this is synchronous as 



1 ] df 



dr 



J 2 2M ^ 2 
dr^ + — ^ dp . 



r y 1-^"' "■ ' r(r,p) 

So the Schwarzschild metric in synchronous coordinates is 

2M 



(4.119) 
(4.120) 

(4.121) 



ds^ = -dr' + 



r{T,p) 



dp2 + r{T, pf {de^ + sin2(e)d0') , (4.122) 



where r{T^ p) is defined by equation fl4.120p . These coordinates are known as 
Leniaitre or Novikov coordinates [^ Wf\ l52] . 

Finally, before taking the limit, shift coordinates 



e' = 6 



vr 



(4.123) 



so that 



ds' = -dt' + 



2M 
r{T,p) 



dp" + r{T, pf (d0'2 + cos2(e') d(/)2) . (4.124) 



Now to take the limit we take p — )■ 0, ^' — )■ 0, such that 

'2 , 



r{r,p) -^r{T) 



r 



(2M); 



and 
Resulting in 

or explicitly, 

ds^ = - 



cosf^') -^ 1. 



2M 



ds^ = -dr^ + ^ dp2 + r(r)2 (d^'^ + d^^) 
r[T) ^ ' 



(4.125) 
(4.126) 

(4.127) 



2 4 

2 . / ^]±\ ' d/ + (2M)i ("^r] ' (d^'2 + d<f) . 



dT^ + 



(4.128) 
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Now in this form we can change the coordinates again to put into a more 
famihar form by rescahng the spatial coordinates: 

ds^ = -dr^ + r-§dp2 + ri (d^'^ + d^') , (4.129) 

which is the Kasner solution with exponents (— |, |, |). 

Equation fl4.129p is a type of Kasner solution (equation (I3.43P ). and this 
particular Kasner solution is also known as the AH I metric [2^ [52] . and 
was initially investigated as a variation on the Schwarzschild metric: 

Then the metric 04.129P can be obtained by performing the coordinate trans- 
formation 



3 \3 /3M 



^=(^) /^; r=[ — ]r-K (4.131) 



3 



2 



4.10. Discussion 



We have developed an ultra-local limit that can be applied to any timelike 
geodesic of an arbitrary spacetime. This construction seems to be the closest 
timelike analogue to the lightlike Penrose limit. Although the construction of 
this ultra-local limit is similar to the Penrose limit, it does not share all the 
familiar properties of that limit, and in particular the hereditary property 
does not follow. 

The result of this limit is a general Bianchi type I metric, which has been 
studied thoroughly in cosmology, both as a simple cosmological model and 
as an ultra-local approximation to more general scenarios. In the context of 
this limit, in contrast to these cosmological settings, there is no clear reason 
to restrict to diagonal metrics, and we have outlined clearly when this is 
a reasonable assumption. Keeping this caveat in mind, the vast amount of 
literature regarding Bianchi type I spacetimes can be a resource of knowledge 
to apply to the ultra-local limit developed herein. 

Working in this general framework we have explored the classical prop- 
erties of the Bianchi type I spacetime, including an examination of Einstein 
equations, covariant conservation, and the Raychaudhuri equation. Further- 
more we have worked through a few simple examples of this limit . 



r CHAPTER 5 

Conclusions 



This thesis has dealt with two research problems in classical general relativ- 
ity, both of which have employed and extended our knowledge of simple, but 
physically useful, exact solutions. These two research problems can be re- 
garded independently, and as such the conclusions at the end of each chapter 
serve to summarise the results for each problem. However, the overall key 
conclusions and some possible extensions are stated below. 



For the first research problem, strong-field gravitational waves were inves- 
tigated. Ways of representing polarisation modes for these waves were ex- 
amined, particularly the previously unaddressed issue of how to write down 
general polarisation modes for the Rosen form. Two forms are commonly 
used for gravitational waves. In the Brinkmann form it is entirely trivial to 
construct any polarisation, whereas it was previously unclear how to do so 
for the Rosen form. 

It was discovered that it was possible to decouple the vacuum equations 
for the Rosen form into parts depending only on an arbitrary walk in some 
appropriately defined polarisation space, and a further differential equation 
relating this walk to the "envelope" function. In (3+1) dimensions this polari- 
sation space is the hyperbolic plane. From this it was possible to write down a 
simple algorithm for constructing polarisation modes. This was furthermore 
extended to arbitrary dimensions. The particularly important, previously 
unknown, circular polarisation was explicitly written down. 

This investigation essentially closes a long-standing gap in our under- 
standing of the Rosen form and strong-field polarisation. An addition that 
could be made is the construction of further specific polarisations, including 
an explicit construction of the elliptical polarisation. 



The second research problem was somewhat more open-ended. An ultra- 
local limit was developed as a close timelike analogue the well-known and 
well studied lightlike Penrose limit. Although this limit seems to be the 
closest timelike analogue to the Penrose limit, this limit does not share all 



69 



70 



the familiar properties of lightlike Penrose limits. 

The result of this ultra-local limit is a general Bianchi type I spacetime. 
Such spacetimes have been most frequently been studied both as simple cos- 
mological models and in approximations to more general cosmological sce- 
narios, where gradient terms may become negligible. In these cosmological 
contexts it is most common to examine diagonal metrics with vacuum, dust, 
or perfect fluid scenarios. Such scenarios do not cover the full range of gen- 
erality desired for our limit, as we saw from an examination of the effect the 
limit had on various scenarios for matter content. Thus we have explicitly set 
forward the necessary and sufficient conditions under which a Bianchi type 
I spacetime is diagonal. 

We have examined several aspects of the Bianchi type I spacetime keep- 
ing full generality. This includes various efforts to write down the useful 
geometrical tensors for the Bianchi type I spacetime as well as an exami- 
nation and interpretation of the Einstein equations, covariant conservation, 
and Raychaudhuri equation. Furthermore, we have looked at a few simple 
examples of how this limit is approached in various scenarios. 

The work on this problem does raise some further avenues for investi- 
gation. Penrose limits have generated considerable interest, and have been 
very well studied in a number of speciflc spacetimes. We have already worked 
through a few simple examples of this limit, and this work could probably 
be usefully extended, hopefully flnding a small number of limits for various 
physically relevant spacetimes. The question still remains about exactly how 
much useful information is retained by this ultra-local limit, a question that 
has been extensively addressed in the case Penrose limits over many years. 
One goal here would be to completely characterise at least some aspects of 
the original spacetimes by properties of the various Bianchi type I spacetimes 
that result from the limits of such spacetimes. 



nAPPENDIX A- 



Publications 



A.l. General polarization modes for the Rosen gravitational wave 



Bethan Cropp and Matt Visser 

Strong-field gravitational plane waves are often represented in either the 
Rosen or Brinkmann forms. While these two metric ansatze are related 
by a coordinate transformation, so that they should describe essentially the 
same physics, they rather puzzlingly seem to treat polarization states quite 
differently. Both ansatze deal equally well with + and x linear polarizations, 
but there is a qualitative difference in the way they deal with circular, ellip- 
tic, and more general polarization states. In this article we will develop a 
general formalism for dealing with arbitrary polarization states in the Rosen 
form of the gravitational wave metric, representing an arbitrary polarization 
by a trajectory in a suitably defined two dimensional hyperbolic plane. 

Published: Classical and Quantum Gravity 27 (2010) 165022. 
arXiv: gr-qc/1004.2734 
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A. 2. Polarization modes for strong-field gravitational waves 



Bethan Cropp and Matt Visser 

Strong-field gravitational plane waves are often represented in either the 
Rosen or Brinkmann forms. These forms are related by a coordinate trans- 
formation, so they should describe essentially the same physics, but the two 
forms treat polarization states quite differently. Both deal well with linear 
polarizations, but there is a qualitative difference in the way they deal with 
circular, elliptic, and more general polarization states. In this article we will 
describe a general algorithm for constructing arbitrary polarization states in 
the Rosen form. 

To be published in the Journal of Physics: Conference Series. 

Proceedings of the Spanish Relativity Meeting (ERE2010). 

arXiv: gr-qc/1011.5904 
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A. 3. Any spacetime has a Bianchi type I spacetime as a limit 



Bethan Cropp and Matt Visser 

Pick an arbitrary timelike geodesic in an arbitrary spacetime. We demon- 
strate that there is a particular limiting process, an "ultra-local limit", in 
which the immediate neighbourhood of the timelike geodesic can be "blown 
up" to yield a general (typically non-diagonal) Bianchi type I spacetime. This 
process shares some (but definitely not all) of the features of the Penrose 
limit, whereby the immediate neighbourhood of an arbitrary null geodesic is 
"blown up" to yield a pp-wave as a limit. 

Published: Classical and Quantum Gravity 28 (2011) 055007. 
arXiv: gr-qc/1008.4639 
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